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AXISYMMETRIC  FILAMENTARY  STRUCTURES 


By  A.F.  Fraser,  P.R.  Preiswerk, 
M.D.  Benton,  and  O.R.  Burggraf 
Astro  Research  Corporation 


SUMMARY 

The  theory  of  filamentary  axisymmetric  structures  is  broad¬ 
ened  to  include  surface  loads  other  than  normal  pressure. 
Structures  with  two  sets  of  symmetrically  disposed  fiber  are 
considered  in  detail,  and  force  transfer  between  filaments  is 
accounted  for  in  the  theory. 

The  governing  equations  are  derived,  and  isotensoid  surface 
shapes  are  determined  and  classified  for  a  wide  range  of  load 
conditions.  The  governing  equations  are  also  applied  to  pro¬ 
blems  where  the  surface  shape  and  loading  are  prescribed  and  the 
filament  geometry  and  load  variation  in  the  surface  are  to  be 
determined. 

The  general  nature  of  filamentary  structures  is  discussed, 
in  depth,  and  useful  qualitative  results  are  obtained  for  multiple 
layer  nets. 

A  non-axi symmetric  isotensoid  (or  isocompressoid)  is  obtained 
as  a  special  limiting  case.  The  shape  is  a  cylindrical  parabo¬ 
loid,  and  the  load  is  equivalent  to  a  gravity  force  field. 
Application  in  civil  engineering  structures  is  indicated. 
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INTRODUCTION 


The  observed  strength  of  many  materials  is  often  greatly 
increased  if  they  are  tested  in  the  form  of  thin  filaments.  For 
example,  continuous  glass  filaments  can  be  produced  with  repro¬ 
ducible  tensile  strengths  ranging  from  200  000  to  700  000  psi. 

Bulk  glass,  on  the  other  hand,  is  usually  limited  to  tensile 
strengths  in  the  order  of  10  000  to  50  000  psi. 

The  high  strength  of  filaments  is  being  exploited  at  an 
increasing  rate  in  the  aircraft,  space,  and  missile  industries. 

The  volume  of  recently  published  information  relating  to  filament 
materials  is  a  good  indication  of  the  current  level  of  effort  ex¬ 
pended  in  both  industrial  and  Government  laboratories  determining 
material  properties  and  performing  design  and  stress  analyses. 

In  the  analysis  and  design  areas,  there  are  two  differeht  but 
overlapping  approaches;  one  is  associated  with  composite-material 
structures  and  the  other  with  filament-dominated  structures.  By 
way  of  definition,  in  the  former,  the  filaments  are  embedded  in 
a  structural  matrix  which  to  a  significant  extent  affects  the 
performance  of  the  material  in  structural  applications.  In  the 
latter,  whether  or  not  there  is  a  matrix,  the  stiffness  and 
strength  of  the  filaments  completely  dominate  the  structural  per¬ 
formance  of  the  material. 

Methods  of  analysis  of  composite  materials  are  based  on 
anisotropic  elasticity  or  shell  theories  (see  e.g.  Ref.  1) .  The 
accuracy  of  these  analyses  depends  on  the  care  taken  in  obtaining 
the  moduli  and  strength  of  the  composites.  A  wealth  of  published 
information  is  available  in  this  field  (see  e.g..  Ref.  2).  The 
accuracy  of  the  analysis  of  filament-dominated  structures  depends 
on  the  degree  to  which  the  filaments  actually  dominate  the  re¬ 
sponse  of  the  structures  analyzed. 

As  methods  of  obtaining  composite  material  properties  im¬ 
prove  and  as  techniques  of  including  effects  other  than  axial 
filament  stiffness  and  strength  enter  into  filament-dominant  anal¬ 
yses,  the  two  methods  should  coalesce  to  the  point  where  materials 
that  are  marginally  in  one  class  or  another  may  be  treated  by 
either  method.  The  work  reported  in  this  paper  broadens  the  anal¬ 
ysis  of  structures  made  of  filament-dominated  materials  to  include 
the  effect  of  force  transfer  between  filaments  or  interfilament 
forces . 


In  order  to  identify  the  value  of  accounting  for  interfila 
ment  forces,  it  is  worthwhile  to  trace  the  recent  development  of 
the  analysis  of  structures  in  which  filament  strength  and  stiff¬ 
ness  play  dominent  roles.  This  research  is  almost  exclusively 
confined  to  axisymmetric  membrane  or  net  structures  which  are 
subjected  to  axisymmetric  surface  and  edge  loads. 

Examination  of  the  equation  for  membrane  stresses  in  shells 
of  revolution  subject  to  uniform  internal  pressure,  p  ,  revea  s 
one  of  the  major  early  stimulants  to  filamentary-structure  de¬ 
sign.  For  this  loading  (see  Ref.  3)  . 


where  Nn  is  the  circumferential  load  per  inch  m  the  membrane, 
and  ri  and  r2  are  the  principal  radii  of  curvature  of  the 
surface,  as  defined  in  Reference  3.  It  is  apparent  that  for 


)  • -  -  _ 

r  /r  =  2  throughout  the  shell,  there  is  no  stress  in  the  cir¬ 
cumferential  direction.  This  led  to  the  supposition  that  the 
phenomenal  strength  of  thin  filamentary  materials  could  be  used 
in  engineering  applications.  For  the  surface  defined  by  r^/r 
=  2  filaments  positioned  on  meridional  lines  would  carry  the 
meridional  stresses  with  uniform  filament  tension  and  no  load 
carrying  capacity  in  the  circumferential-  direction  would  be  re¬ 
quired  except  at  the  edges  of  the  surface. 

This  concept  was  applied  in  the  design  of  parachutes  as 
early  as  1923  (Ref.  4),  and  more  recently  m  the  design  develop¬ 
ment  of  radial  cord  tires.  In  1959,  end  closures  for  filament- 
wound  pressure  vessels  (rocket  chambers)  were  designed  (Refs.  5, 
14  and  15)  In  1960  a  "radial  cord"  toroidal ,  deployable  space- 
stitton  model  was  designed  (Ref.  16)  .  and  48-inch  and  24-foot 
diameter  engineering  test  models  were  constructed  and  subsequent¬ 
ly  evaluated  (Refs.  17,  18,  and  19).  In  1965  an  airlock  design 
for  space  application  was  proposed  (Ref.  6).  In  all  these  app  l 
cations  a  membrane  is  required  to  contain  internal  pressure  and 
to  introduce  the  loads  into  the  filaments. 

With  the  exception  of  the  airlock  and  the  toroidal  configur 
ations,  the  disadvantages  of  filament  buildup  at  the  apex  o  e 
membrane  led  to  examination  of  other  surface  shapes  in  *ni 

form  tension  in  filaments  placed  on  non-meridional  lines  would 
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contain  the  internal  pressure.  The  problem  was  formulated  in 
Reference  5  for  the  simple  case  of  two  sets  of  filaments  sym¬ 
metrically  disposed  with  respect  to  meridional  lines,  and  the  re¬ 
sulting  shapes  were  determined  in  Reference  7  for  various  edge 
conditions . 

The  isotensoid  concept,  that  of  placing  filaments  on  sur¬ 
faces  of  revolution  so  that  design  loads  are  precisely  balanced 
by  axial  and  uniform  stresses  in  the  filaments,  was  then  put  on 
a  firm  basis  in  1962  with  the  appearance  of  Reference  8.  In  this 
reference,  the  complete  theory  of  pressurized  filament -wound 
shells  was  formulated;  a  systematic  study  of  isotensoid  shapes 
was  made,  and  the  basis  for  analysis  which  includes  centrifugal 
loads  caused  by  rotation  was  explored.  Also,  the  equations  were 
developed  for  nonuniform  pressure  and  multiple-layer  shells,  and 
a  wide  range  of  isotensoid  structures  was  described  in  Reference 
20.  The  isotensoid  shapes  determined  in  Reference  8  made  pos¬ 
sible  the  solution  of  the  indirect  design  problem  in  some  cases. 
That  is,  given  a  mission  objective,  if  one  of  the  shapes  result¬ 
ing  from  the  analysis  would  perform  the  mission,  the  problem  was 
solved. 

The  following  year,  the  problem  of  designing  pressurized 
rotating  net  structures  was  solved  and  published  in  Reference  9. 
The  problem  was  completely  defined  and  all  possible  shapes 
determined  and  reported.  It  was  shown  that  filament  tension  of 
necessity  varies  in  a  spinning  pressurized  net  provided  that 
there  is  no  force  transfer  between  filaments  at  crossover  points. 

Concurrent  with  this  work,  another  approach  to  the  design  of 
pressure  containers  was  taken  and  reported  in  Reference  10. 

Here,  the  direct  problem  was  solved  by  allowing  a  continuous  dis¬ 
tribution  of  filament  angles  or  layers  of  filaments.  Using  the 
results  of  Reference  10,  one  can,  in  principle,  specify  membrane 
geometry,  pressure  loading,  and  isotensoid  filament  stresses  and 
determine  the  required  distribution  of  filaments  in  the  membrane. 
Practical  design  considerations  dictate  the  continuous  distribu¬ 
tion  be  replaced  by  a  finite  set,  thus  approximations  are 
required. 

In  summary,  the  principal  results  of  published  investiga¬ 
tions  into  filament-dominated  designs  can  be  found  in  References 
8,  9,  and  10.  Reference  8  identifies  the  isotensoid  shapes  pos¬ 
sible  for  two-family  symmetrically  wrapped  pressurized  shells  of 
revolution  and  presents  solutions  for  single-layer  (two-family) 
and  multiple-layer  (more  than  two-family)  isotensoids.  Reference 
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9  identifies  shapes  for  spinning  pressurized  nets  formed  by  fila¬ 
ments  symmetrically  disposed  on  a  shell  of _ revolution  and  pre- 
sents  solutions  for  these  shapes  and  equations  for  filament  ^ 
tension  variation.  Reference  10  gives  a  method  of  solution  and 
equations  defining  the  filament  wrap  pattern  for  isotensoid  de¬ 
signs  of  pressurized  shells  of  revolution  when  the  shell  geometry 

is  specified. 

Throughout  the  research  reviewed  above,  an  artificial  design 
constraint  has  been  imposed;  that  is,  force  transfer  between 
filaments  at  filament  crossover  points  is  neglected.  In  nets 
this  transfer  is  manifested  in  knot  forces,  and  in  composites  it 
is  manifested  in  self-equilibrated  matrix  stresses.  In  both 
cases  the  interfiber  force  has  the  characteristics  of  shear  loads, 
and  it  has  been  termed  interfiber  shear  or  interlaminar  shear,  etc 
Unfortunately,  this  terminology  carries  the  connotation  of  shear 
in  continuum  mechanics  and  the  fundamental  nature  of  the  force 
transfer  between  filaments  is  obscured.  For  this  reason,  the 
term  » interfilament  force  transfer"  will  be  used  throughout  this 

report. 

The  design  freedom  obtained  by  employing  the  interfilament 
force  transfer  is  exemplified  in  Reference  11  where  the  solution 
is  obtained  for  an  isotensoid  spinning  filamentary  disk.  If  the 
interfilament  force-transfer  capability  is  ignored  m  a  spinning 
disk,  it  becomes  a  degenerate  case  of  the  general  theory  of  Re 
erence  9  for  which  isotensoid  design  is  impossible.  Therefore, 
consideration  of  interfilament  force  transfer  allows  isotensoid 
(and  minimum  weight)  design  where  it  was  not  formerly  possi  e. 
Also,  as  will  be  indicated  in  the  section  entitled  "Fundamental 
Nature  of  Filamentary  Structures",  the  use  of  interfilament  force 
transfer  in  the  design  process  allows  solution  of  problems  where 
the  shape  is  prescribed  without  resorting  to  continuous  distribu¬ 
tions  of  filaments  as  was  necessary  in  the  approach  taken  in  Ref 
erence  10. 

The  work  reported  in  this  paper  introduces  the  interf ilament 
force  transfer  developed  in  Reference  11  into  the  general  ap¬ 
proach  of  References  8  and  9.  This  is  done  for  arbitrary  sym¬ 
metric  loading  of  axisymmetric  nets  formed  by  two  families  (or 
one  layer)  of  filaments  as  defined  in  the  section  entitled 
"Fundamental  Nature  of  Filamentary  Structures".  It  is  shown  that 
isotensoid  solutions  are  possible  for  many  problems  in  the  pres 
ence  of- any  axi symmetric-surface  load-vector  field  provided  the 
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vector  field  is  confined  to  meridional  planes.  Thus  pressure, 
centrifugal  force,  and  axial  surface  load  can  be  accommodated 
with  isotensoid  minimum-weight  design. 

In  the  section  entitled  "Fundamental  Nature  of  Filamentary 
Structures",  the  fundamental  nature  of  string  structures  is  dis¬ 
cussed  qualitatively.  The  governing  equations  for  two-family 
nets  are  derived  in  the  section  entitled  "Analysis  of  Single  Layer 
Surfaces".  These  equations  are  solved  for  a  number  of  load  and 
design  conditions,  and  the  results  are  reported  and  discussed  as 
they  are  obtained.  The  problem  of  determining  the  filament 
orientation  in  the  surface  and  tension  variation  along  filaments 
when  the  surface  shape  is  prescribed,  is  solved.  Four  illustra¬ 
tive  examples  are  given.  Effort  is  then  concentrated  on  isoten¬ 
soid  design,  and  the  equations  are  specialized  to  apply  to  the 
case  of  uniform  filament  tension.  A  wide  range  of  load  conditions 
is  considered  and  methods  of  solution,  classification  of  isoten¬ 
soid  surface  shapes,  and  numerous  specific  examples  are  given.  A 
summary  of  the  results  is  given  in  the  section  entitled  "Closure", 
and  possible  extensions  are  discussed. 
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SYMBOLS 


a  uniform  acceleration  in  negative  z  direction 

b  parameter  defined  by  Equation  (87) 

d  parameter  defined  by  Equation  (86) 

argument  in  elliptic  integrals  defined  in  Equations  (73) , 
1  (77) ,  (81) ,  and  (85) 

l  distance  measured  along  filament 

fit  distance  measured  along  a  filament  between  two  adjacent 

filament  intersections 

m1  mass  per  unit  length  of  filament 

n  total  number  of  filaments  in  a  layer 

p  pressure 

pQ  pressure  at  r  =  rQ 

r  radial  coordinate 

rl'  r2  Princ:>-Pal  radii  of  curvature  of  surface  (circumferential 
and  meridional  respectively) 

rQ  reference  radius 

rffl  parameter  used  in  allowing  radial  coordinate  to  approach 

infinity 

y  resulting  cartesian  coordinate  when  r^,  — ►  °°  (r  =  r^  +  y) 

z  axial  coordinate 

A  axial  acceleration  load  parameter  (Equation  12) 

A,B,C,D  solution  types  defined  by  Equation  (67) 

A  A/r0 
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constants 


D  denominator  function  (Equation  88) 

E  ,  E  (kjj  elliptic  integrals  of  the  second  kind 


Fn'Fm'Fc 


vector  force  field  acting  on  a  filament  (lb/in) 
force  at  filament  intersections  (lb) 

components  of  F  in  normal,  meridional  and  circumferen¬ 
tial  directions 


F(f^,ki)  elliptic  integral  of  the  first  kind 
G  parameter  defined  by  Equation  68 

K  pressure  load  parameter  (Equation  12) 

K(k^)  complete  elliptic  integral  of  the  first  kind 

L  vector  surface  load  (lb/in  ) 

Lr|/Lrn  components  of  L  in  normal  and  meridional  directions 

^  (lb/in2) 


N 

Q 

R 

T 


T 


o 


X 


X2'X3 


Y  (X) 


Z 


numerator  function  (Equation  88) 
function  defined  by  Equation  (22) 

R  =  r/rQ 

filament  tension 

filament  tension  at  r  =  rQ 

X  =  R2  =  (r/rQ)2 

roots  of  Equation  (66)  given  in  Equation  (67) 
function  defined  in  Equation  (66) 
z  =  z/rc 
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angle  between  meridional  tangent  and  z  axis 

value  of  a  at  r  =  rQ 

angle  between  filament  and  meridian 

value  of  (3  at  r  =  rQ 

longitudinal  coordinate 

argument  in  elliptic  integrals  defined  in  Equations  (72) , 
(76) ,  (80) ,  and  (84) . 

angular  velocity 

spin  load  parameter  (Equation  12) 


FUNDAMENTAL  NATURE  OF  FILAMENTARY  STRUCTURES 


In  order  to  discuss  the  fundamental  nature  of  filamentary 
structures,  it  is  necessary  to  develop  clear  and  concise  defini¬ 
tions.  In  what  follows,  all  discussion  is  restricted  to  axisym- 
metric  surfaces  and  axisymmetric-applied  surface  and  edge  loads 
both  of  which  have  no  components  normal  to  the  meridional  planes. 

The  equation  of  the  surface  is  expressed  in  the  cylindrical 
polar  coordinates  of  Figure  1  as 


z  =  z  (r) 
or  Z  =  Z  (R) 

where  Z  =  z/rQ  and  R  =  r/rQ  ,  and  rQ  is  some  reference  radius. 
The  external  load  vector  field  shown  in  Figure  1 

L  =  L(r) 


has  no  component  normal  to  meridional  planes  and  is  independent  of 
the  longitude ,  I  . 

With  reference  to  Figure  2,  a  filament  set  or  family  is  a 
collection  of  filaments  which  on  parallel  circles  form  equal 
angles,  p  ,  with  meridional  lines  on  the  surface.  Within  a  fila¬ 
ment  set  p  =  P(r)  as  shown  in  Figure  2a.  To  each  filament  set  or 
family  at  P  ,  there  corresponds  a  complementary  set  at  -p  as 
shown  in  Figure  2b.  These  two  sets  together  are  called  a  layer 
(Fig.  2c) .  It  is  apparent  from  the  symmetry  of  the  system  under 
consideration  that  filament  sets  always  appear  in  complementary 
pairs  forming  a  layer.  In  some  cases  a  single  filament  will  be¬ 
long  to  both  sets  in  a  layer,  as  shown  in  Figure  2d.  It  is 
apparent  that  filaments  in  a  set  do  not  intersect  except  on  unique 
parallel  circles  where  p  =  tt/2  ,  and  on  these  circles,  each  fila¬ 
ment  is  leaving  one  set  of  a  layer  and  entering  the  complementary 
set.  Filament  sets  at  p  =  tt/2  and  p  =  0  are  considered  to  be 
degenerate  layers  in  which  complementary  sets  each  include  one 
half  the  filaments  in  the  layer,  and  both  sets  are  the  same. 
Filamentary  shells  or  nets  can  be  designed  with  one  or  more 
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layers;  the  single  layer  nets  will  be  discussed  first. 

Single  Layer  Nets.  -  Consider  a  single  layer  net  as  shown  in 
Figure  2c  and  for  the  moment  let  the  surface  load  vector  be 
directed  normal  to  the  surface  formed  by  the  net.  Further,  assume 
that  the  filaments  in  a  set  are  free  to  slide  relative  to  the 
other  set  in  the  layer.  The  net  will  form  an  axisymmetric  surface 
shape  which  depends  on  the  pressure,  the  inner  and  outer  radii 
(assuming  they  are  fixed) ,  the  length  of  the  filaments,  and  the 
total  angle,  §  , traversed  by  each  filament  as  it  is  traced  from  the 
inner  radius  to  the  outer  radius.  The  filaments  will  be  stressed 
uniformly  and  will  trace  out  geodesic  curves  on  the  surface  form¬ 
ed.  These  facts  are  established  in  Reference  8.  The  net  con¬ 
sidered  is  clearly  a  mechanism  which  assumes  different  isotensoid 
shapes  which  will  depend  on  how  the  surface  load  vector  varies. 

Now  relax  the  condition  that  the  surface  load  vector  be  nor¬ 
mal  to  the  surface  and  allow,  for  example,  centrifugal  surface 
forces.  (This  simple  terminology  is  used  in  place  of  the  more 
exact  expression  "centrifugal  body  force  per  unit  area  of  sur¬ 
face".)  The  surface  shape  and  the  filament  paths  will  change.  In 
fact,  the  filament  geometry  in  the  new  surface  will  be  such  that 
the  surface-load-vector  field  can  be  projected  at  each  point  into 
two  components;  one  in  the  direction  of  the  principal  normal  of 
the  filament  and  the  other  in  the  direction  of  the  tangent  to  the 
filament.  The  latter  surface  force  component  causes  the  tension 
in  each  filament  to  vary  along  its  length.  These  facts  are  estab¬ 
lished  in  Reference  9  and  account  for  the  fact  that  isotensoid  de¬ 
sign  is  not  possible  under  the  load  and  net  conditions  stipulated. 

Now  assume  that  at  each  junction  among  filaments  in  separate 
sets  in  the  layer  there  is  either  a  knot  or  a  dab  of  glue.  As  a 
consequence,  if  there  is  any  tendency  of  the  filaments  in  comple¬ 
mentary  sets  to  slide  relative  to  one  another,  a  constraining 
force  will  result  in  the  knot  or  glue  and  the  sliding  motion  will 
be  arrested.  The  axisymmetry  of  the  net  and  the  loading  do  not 
allow  the  force  transferred  from  one  filament  to  the  other  to 
have  a  component  in  the  meridional  direction  of  the  surface.  The 
interfilament  force  transfer  between  filaments  in  complementary 
sets  of  a  layer  must,  therefore,  be  in  the  circumferential  direc¬ 
tion  as  shown  in  Figure  3.  It  is  termed  Fk  .  Thus,  if  there  is 
any  tendency  for  slip  to  occur  between  filament  sets  in  a  layer, 
even  though  the  applied  load  is  restricted  to  lie  in  meridional 
planes,  each  filament  will  experience  a  load  vector  with  a 
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component,  F^  ,  out  of  the  meridional  plane  at  each  of  the  inter¬ 
sections  with  fibers  from  the  complementary  set.  These  loads, 

Fk  ,  are  internal  and  self  equilibrating. 

Now  consider  the  effect  of  interfilament  force  transfer  on 
the  two  problems  discussed  above.  For  the  case  of  surface  load¬ 
ing  normal  to  the  surface,  imagine  that  a  pressure  has  been 
applied  to  a  system  in  which  sliding  can  occur  and  the  appropriate 
isotensoid  shape  has  been  assumed.  Now  imagine  that  the  filaments 
are  either  knotted  or  glued  at  each  intersection.  It  is  clear 
that  the  surface,  the  filament  geometry  in  the  surface,  and  the 
stress  level  in  the  filaments  is  undisturbed.  Consequently,  the 
presence  of  knots  does  not  increase  the  number  of  design  parame¬ 
ters  for  the  geometry  and  load  considered. 

In  the  second  case  where  the  surface-force-vector  field  is 
not  normal  to  the  surface,  filament  force  transfer  can  be  used  to 
much  greater  advantage.  It  has  been  established  (Reference  8) 
that  isotensoid  design  is  not  possible  in  the  absence  of  inter- 
filament  force  transfer  when  the  surface  load  vector  is  not  nor¬ 
mal  to  the  surface.  Consider  now  a  design  where  interfilament 
force  transfer  exists  and  non-normal  surface  load  is  applied. 

The  possibility  exists  that  the  filament  geometry  (and  therefore 
the  filament  intersections)  and  the  surface  shape  can  be  deter¬ 
mined  so  that  the  vector  sum  of  external  load  per  filament  and  the 
interfilament  load  has  its  direction  along  ;the  principal  normal  to 
the  filament  path  in  the  surface.  The  resulting  net  configuration 
will  be  isotensoid  under  the  applied  loads.  Then  the  use  of  fila¬ 
ment  force  transfer  in  the  design  stage  will  allow  optimum  design 
and  minimum  weight  in  single  layer  nets  when  the  surface  load  vec¬ 
tor  field  is  not  normal  to  the  surface.  The  examination  of  this 
problem  is  the  subject  of  a  large  part  of  the  remainder  of  this 
report. 

Multiple  Layer  Nets.  -  Consider  a  system  of  two  layers  form¬ 
ing  an  axisymmetric  surface  and  subjected  to  the  type  of  general 
symmetric  loads  defined  above.  Forces  tangent  to  the  surface  (in 
addition  to  bearing  forces)  can  be  transferred  between  layers, 
and  in  addition  to  the  transfer  of  a  part  of  the  tangential  com¬ 
ponent  of  the  surface  load  to  the  inner  layer,  a  self —equilibrated 
system  of  loads  between  layers  similar  to  that  between  sets  may 
arise.  Symmetry  considerations  restrict  the  self-equilibrated 
load  system  between  layers  to  be  directed  along  meridians  of  the 
surface.  The  load  on  a  filament  is  then  the  vector  sum  of  the 
applied  load,  a  circumferential  load  due  to  the  self-equilibrated 
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forces  between  complementary  sets,  and  a  meridional  load  due  to 
the  self-equilibrated  forces  between  layers.  The  possibility 
then  exists  that  the  shape  and  surface  load  be  specified  and  the 
filament  paths  in  both  layers  determined  so  that  the  resultant 
filament  loads  are  directed  along  the  principal  normals  of  the 
filaments.  Under  these  conditions,  isotensoid  design  is  possible 
with  a  finite  number  of  layers  even  when  the  geometry  is  speci¬ 
fied.  In  the  absence  of  interfilament  force  transfer,  this  pro¬ 
blem  has  been  solved  only  for  normal  surface  loads,  and  an 
infinite  number  of  layers  was  required  (Ref.  10) . 

A  basic  difference  between  single-layer  nets  and  multiple- 
layer  nets  has  implicitly  manifested  itself  in  this  section. 

That  is,  single-layer  nets,  even  with  pinned  intersections,  are 
mechanisms  and  the  shape  taken  by  them  is  strongly  load  dependent. 
On  the  contrary,  multiple-layer  nets,  if  knotted  or  glued,  are 
redundant  structures,  and  fundamentally  different  design  freedom 
can  be  achieved. 
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ANALYSIS  OF  SINGLE  LAYER  NETS 


Consider  equilibrium  of  a  net  formed  by  one  layer  (Fig.  2c) 
of  complementary  sets  as  defined  in  the  section  entitled  "Funda¬ 
mental  Nature  of  Filamentary  Structures".  The  filaments  making 
up  the  sets  are  flexible  but  essentially  inextensible.  Each 
intersection  between  filaments  in  complementary  sets  is  pinned 
either  by  a  knot  or  a  dab  of  glue.  The  surface  load  is  a  vector 
field  with  components  confined  to  meridional  planes  (Fig.  1)  and 
is  independent  of  longitudinal  position  on  the  surface. 

Let  F  be  the  load  vector  per  unit  length  acting  on  a  fila¬ 
ment.  The  components  of  F  in  the  normal  and  meridional  direc¬ 
tions  of  the  surface  are  related  to  the  similarly  directed 
component  of  surface  loading,  L  ,by 


_  2TTr  cos  ft 


n 


Li 


(1) 


where  ft  is  the  filament  angle  defined  in  Figure  4,  n  is  the 
total  number  of  filaments  in  the  surface,  and  the  subscript,  i  , 
denotes  either  the  meridional  or  normal  direction. 

In  order  to  obtain  the  component  of  F  in  the  longitudinal 
direction,  assume  that  filament  crossover  points  or  intersections 
are  sufficiently  close  that  the  filament  transfer  force  can  be 
averaged  between  intersections.  Then  an  interfilament  force  per 
unit  length,  Fc  ,  can  be  defined  where  Fc  At  =  F^  and  At  is 
the  distance  along  a  filament  between  intersections.  The  total 
load  per  unit  length  of  filament  is  the  vector  sum  of  the  ortho¬ 
gonal  components. 

Equilibrium  of  a  surface  element  is  enforced  by  considering 
equilibrium  of  filaments  in  one  set  and  accounting  for  the  pre¬ 
sence  of  the  other  set  through  the  filament  transfer  force,  Fc  . 
To  this  end,  consider  equilibrium  of  a  filament  of  length,  dl  , 
(Fig.  4)  where  along  a  filament,  the  following  differential  rela¬ 
tions  hold 
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dr  =  d-t  cosg  sina 


(2) 


d$  =  d-t  r--'-  sing 
dz  -  -d-t  cos  3  cosa 


where  a  is  the  angle  between  a  meridian  on  the  surface  and  the 
z-axis  as  shown  in  Figure  4. 

With  T  =  T(r)  the  tensile  force  in  the  filament,  three 
independent  equilibrium  equations  can  be  written;  equilibrium  of 
forces  parallel  to  the  filament,  equilibrium  of  forces  in  the  z- 
direction,  and  equilibrium  of  torques  about  the  z-axis.  These 
equations  read  respectively: 

dT  -  F  sing  d-t  -  F  cosg  d-t  =  0  \ 

c  m  I 

d  (T  cosg  cosa)  -  Fn  sina  d-t  -  Fm  cosa  d-t  =  CM  (3) 

d(rT  sing)  -  Fcr  d-t  =  0  I 

Equation  2  can  be  used  to  eliminate  d-t  from  Equation  (3)  . 
With  r  chosen  as  the  independent  variable.  Equations  (3)  then 
read: 


ZI.  -  Fc  tang  csca  -  Fm  csca  =  0  (4) 


d(T  cos g  cosa)  _  Fn  secg  -  F  secg  etna  =  0  (5) 

dr 


d  (rT  sing)  -  f  r  secg  csca  =  0 

-i  c 

dr 


(6) 


This  system  of  first-order  differential  equations  governs 
the  design  and  analysis  of  single-layer  nets.  The  functions, 
g(r)  and  a (r)  ,  can  be  thought  of  as  the  intrinsic  coordinates 
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of  the  filament  paths  on  the  surface  and  the  surface  shape,  re¬ 
spectively.  They  are  related  to  the  more  practical  design  func¬ 
tions,  § (r)  and  z (r)  ,  by  two  first-order  differential 
equations  obtained  by  eliminating  d-i  from  Equation  (2)  . 


=  r--*-  tanp  csca  (7) 

dr 


—  =  -  etna  (8) 

dr 


Equations  (4)  through  (8)  have  been  examined  in  detail  for 
the  special  case  of  combined  pressure  and  centrifugal  loading 
with  no  force  transfer  between  filaments  (Fc  =  0)  in  Reference 
9,  and  for  pressure  loading  (Fm  =  FQ  =  0)  in  Reference  8  where 
an  arbitrary  number  of  layers  was  considered.  In  neither  case 
could  the  shape  be  prescribed.  As  pointed  out  in  the  section 
entitled  "Fundamental  Nature  of  Filamentary  Structures",  force 
transfer  between  filaments  in  separate  sets  of  a  layer  allows 
specification  of  the  surface  shape  provided  filament  tension  may 
vary  along  a  filament.  This  extends  the  design  freedom  in  pro¬ 
blems  of  the  type  treated  in  Reference  9.  Further,  the  inclusion 
of  filament  force  transfer  allows  isotensoid  design  when  the  sur¬ 
face  load  vector  is  not  normal  to  the  surface.  Thus,  the  results 
of  Reference  8  can  be  extended  to  determine  isotensoid  shapes  for 
arbitrary  symmetric  surface  loading. 

Two  distinct  classes  of  problems  are  thus  identified  for 
single  layer  nets : 

(1)  The  shape  and  surface  load  are  prescribed  and  the  fila¬ 
ment  tension,  T(r)  ,  and  the  filament  geometry,  i (r)  ,  are  to 

be  determined. 

(2)  The  surface  load  and  uniform  filament  tension  are  pre¬ 
scribed  and  the  surface  shape,  z (r)  ,  and  the  filament  geometry, 

$  (r)  ,  are  to  be  determined. 

Both  of  these  problems  are  treated  below.  Under  certain  condi¬ 
tions  a  third  type  of  problem,  not  considered  here,  must  be 
addressed.  If  the  interfilament  forces  become  excessive,  failure 
can  occur  between  sets  in  a  bonded  design  or  in  knots  in  a  knotted 
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design.  Therefore,  in  any  design  application,  an  allowable  must 

be  established  for  F  and  care  taken  that  it  is  not  exceeded. 

c 


Before  solutions  for  the  two  problem  types  can  be  developed, 
the  explicit  dependence  of  Fm  and  Fn  on  the  independent 
variable,  r  ,  must  be  stipulated.  Three  specific  load  cases  are 
considered;  surface  pressure,  p  ,  centrifugal  surface  load 
caused  by  angular  rotation,  uu  ,  about  the  axis  of  symmetry,  and 
a  d'Alembert  surface  load  in  the  positive  g-direction  caused  by  a 
uniform  acceleration,  "a"  ,  in  the  negative  z-direction. 


For  pressure,  p(r) 


Fn  =  cos^ 


F  =  0 
m 


(9) 


For  radial  acceleration  caused  by  spin  velocity,  uu 


2 

F  =  m’o)  r  coscc 
n 

2 

F  =  -  m‘  uo  r  sina 


(10) 


For  axial  acceleration,  "a"  ,  in  the  negative,  z-direction 


Fn  =  m 1  a  sina  j 

(id 

F„,  =  m'a  coscc  I 

m  | 


where  m'  is  the  filament  mass  per  unit  length.  Note  that  the 
presence  of  a  dense  fluid  contained  within  the  surface  and  under¬ 
going  acceleration  and  rotation  may  be  accounted  for  through  the 
pressure,  p(r)  .  It  is  advantageous  at  this  point  to  nondimen- 
sionalize  the  equations.  Let  TQ  and  pQ  be  the  filament 
tension  and  pressure  load  at  some  reference  radius,  rQ  ,  and  let 
the  following  dimensionless  quantities  be  defined 
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t 

t 


Pressure  parameter 


Spin  parameter 


2rrr  z 

K  -  - —  p^ 

nT  ° 


m 1  uo^r  ^ 

fi  =  - 2_ 


Axial  acceleration  parameter  A  = 


m  ar 


T, 


Radial  coordinate 


R  =  r/r, 


k-  (12) 


The  filament  loads  per  unit  length,  Fm  and  Fn  ,  are  convenient¬ 
ly  expressed  by  use  of  the  definitions  (12)  and  summing  Equations 
(9)  through  (11) . 


T 

F  =  _o  (-Q  R  sinOC  +  A  coscx)  \ 

III  y'  I 

°  / 

[  (13) 

m  1 

Fn  =  —  (KR  p/pQ  cosp  +  Q  R  cosa  +  a  sina)  I 

r  / 


Equation  (13)  for  Fm  and  Fn  can  be  used  in  Equations  (4) 
through  (6)  to  obtain  the  nondimensional  governing  equations  for 
single-layer  surfaces  of  revolution  loaded  by  internal  pressure, 
axial  acceleration,  and  spin. 


d(T/To> 

dR 


tan (3  csca  +  Q  R  -  A  etna 


0 


(14) 


d(cosp  cosa  t/T  ) 
dR 


-  KR  p/pQ  -  A  seep  csca  =  0 


(15) 
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nT^TJTA 


(16) 


d  (R  sinf3  T/Tq)  _  /Fg  rQ 
dR  \  Tq 

For  convenience.  Equations  (7)  and  (8) ,  which  are  the  differ¬ 
ential  relations  among  the  intrinsic  coordinates  ( (3  and  a)  , 
and  the  more  practical  design  functions  (z  and  I)  ,  are  re¬ 
corded  in  nondimensional  form.  Let  z/rQ  =  Z  in  Equations  (7) 
and  (8) ,  then 

M.  =  R-^  tanP  csca  (17) 

dR 


=  -  etna  (18) 

dR 


R  secP  csca  =  0 


The  remainder  of  this  work  is  devoted  to  the  solution  of 
these  governing  equations  for  two  types  of  problems: 

(1)  Surface  shape  and  loading  prescribed,  filament  tension 
and  filament  angle  to  be  determined,  and 

(2)  Uniform  tension  and  loading  prescribed  and  surface 
shape  to  be  determined. 

In  problems  of  the  first  type,  it  is  possible  to  apply  edge  bound¬ 
ary  conditions  on  loads  by  use  of  either  Equation  (23)  or  (35) 
below.  However,  it  is  more  convenient  in  the  problems  treated  to 
specify  the  filament  angle  at  the  edge  and  determine  the  edge 
filament  loads  parametrically  in  terms  of  the  specified  angle. 

In  problems  of  the  second  type,  the  shape  and  edge  position  are 
unknown  at  the  outset,  consequently  it  is  convenient  again  to  use 
the  filament  angle  pQ  at  R  =  1  as  a  boundary  condition  and  to 
determine  the  shapes  with  pQ  as  a  parameter.  The  shapes  thus 
determined  may  terminate  at  edges  where  p  ^  tt/2  .  Then,  with 
the  shape  and  filament  angle  known  at  the  edge,  the  edge  load 
conditions  necessary  to  support  the  filament  loads  may  then  be 
determined. 


The  Direct  Problem  -  Shape  Prescribed 
Barring  the  exceptional  case  when  the  prescribed  shape  is 
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an  isotensoid  shape  for  the  loads  applied,  the  tension  in  the 
filaments  of  a  single  layer  surface  will  be  functions  of  r  . 

With  the  loads  (K,  A,  Q)  and  shape  (a)  prescribed,  Equations 
(14),  (15),  and  (16)  must  be  used  to  determine  the  unknown  func¬ 
tions,  F  (r)  ,  p (r)  ,  and  T(r)  .  The  function,  F  (r)  ,  is  easily 

eliminated  from  Equations  (14)  and  (16),  leaving  two  governing 
equations 


sinP 


d (R  sinP  T/Tq) 


dR 


=  R 


d (T/T0) 
dR 


+  fi  R  -  AR  etna 


(19) 


d(cosP  coset  T/T  ) 
_ 


KR  p/Pq  +  A  secP  esca 


(20) 


For  A  ^  0  ,  the  appearance  of  p  in  the  last  term  in  Equa¬ 
tion  (2)  precludes  its  solution  by  simple  integration.  The  two 
cases,  A  =  0  and  A  ^  0  ,  are  therefore  treated  separately. 

Solution  When  A  =  0.  -  For  A  =  0  ,  solutions  when  they 
exist  can  be  expressed  analytically.  Equation  (20)  can  be 
integrated 


cosp  = 


K  seca 


P/P0  RdR 


As  boundary  condition,  set  P  =  P0  and  sinPQ  =  C  at  R  =  1  and 
for  convenience,  define  the  new  independent  variable,  X  =  R^  . 
Then 


cosP  =  seca 


(21) 


The  right  side  of  Equation  (21)  is  a  known  function  of  the  inde¬ 
pendent  variable,  X  =  R^  =  (r/rQ) ^  ,  which  depends  on  the  shape 
and  prescribed  pressure  load.  Let  it  equal  Q (X)  .  Then 
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Q(x) 


seca 


(22) 


1  -  Cz  cosa 


-?/ 


(p/pD)  dx 


X 


and  Equation  (21)  reads 


=  Q  sec3  (23) 

To 

Equation  (23)  can  be  used  to  eliminate  T/Tq  from  Equation  (19) 
when  A  =  0  with  the  result 


sinp  [R  Q  tanp] 
dR 


[Q  seep]  +  n  R2 


If  this  expression  is  multiplied  by  Q  sec (3  and  the  change  to 
independent  variable,  X  ,  completed,  the  following  expression 
results 


d 

dX 


[X  Q2 


tan2  3]  X  —  [Q2 


sec2 j5]  +  Q  X  Q  secP 


(24) 


Use  of  the  identity  tan2  3  =  sec2 3-1  in  (24)  and  recognition 
that 


d_ 

dX 


[X  Q2  sec2  p]  =  X 


d_ 

dX 


[Q- 


sec2  3] 


+  Q  sec 


results  in  the  following  quadratic  equation  for  sec 3 


sec2  3 


-  Q  X  Q  sec  3 


-St* 


Q2]  = 


0 


(25) 


from  which, 


20 


(26) 


seep 


n  X 


±  V 


nz  xz  +  4 


[x 

dx 


Q2] 


The  sign  in  (26)  is  chosen  so  that  seep  >  1  .  With  this  result 
for  the  filament  geometry,  the  tension  is  given  by  Equation  (23) . 


ft  X  + 


ft  X  +  4 


dX 


[X  Q2] 


(2  7.) 


where  the  sign  is  chosen  as  for  seep  .  The  interfi lament  force 
can  now  be  computed  directly  from  Equation  (16) .  The  filament 
path  in  terms  of  the  polar  coordinates  is  determined  by  integra¬ 
tion  of  Equation  (17)  for  p  given  by  Equation  (26) . 

From  these  results  it  is  apparent  that  not  all  the  parameters 
can  be  chosen  independently.  For  example,  if  Equation  (25)  is 
evaluated  at  X  =  1  ,  the  following  condition  results 


ft  =  1 


d_ 

dx 


[X  Q2] 


(X 


1) 


(28) 


To  obtain  Equation  (28),  note  that  Q(l)  was  determined  from 
Equation  (22)  as 


Q  (1)  =  V1  ~  °2  =  cosPo  (29) 


Equation  (28)  may  be  regarded  as  a  relation  between  PQ  and  TQ 
for  geometry,  load,  and  all  other  parameters  fixed.  Thus  with 
either  TQ  or  pQ  specified  at  r  =  rQ  ,  Equations  (26)  and 
(27)  determine  the  unique  solution  for  tension  and  filament  path 
in  the  specified  surface  subjected  to  the  specified  loads. 

The  minimum  (maximum)  radius  to  which  the  solution  extends 
on  the  specified  shape  occurs  at  p  =  tt/2  or  p  =  0  .  For  the 
filament  tension  to  remain  finite,  when  P  =  tt/2  ,  Equation  (23) 
requires  that 
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OlXmax)  '  0 


(30) 


For  this  condition  Equation  (22)  determines  Xmax  (or 


1 

j  (p/P0)  dX 
•^max 


2 

K 


cosa0 


(31) 


Equation  (31)  defines  the  minimum  radius  for  internal  pressure 
(K  >  0)  and  the  maximum  radius  for  external  pressure  (K  <  0)  . 

Finally,  if  the  filaments  become  tangent  to  the  meridians, 
the  solution  ends.  In  this  case  3=0  and  Equation  (25)  re¬ 
quires  that 

dQ  =  _  H  -For  p  =  o  (32) 

dX  2 


Three  examples  are  given  to  illustrate  the  solution  method. 

Example  1  -  Spherical  Pressure  Vessel.  -  Q  =  0  ,  P  =  PQ 
The  shape  of  the  sphere  is  defined  by 


Hence 


and 


cosa 


x. 

X  ^ 


[x  Q2] 

dX 


=  K 


-  £  (1  -  X) 
2 


With  this  result,  the  boundary  condition.  Equation  (28) ,  gives 
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K 


1 


=  - .  :  =  =  seep 

V 1  -  C2  ° 

which  with  the  definition  of  K  in  Equation  (12)  determines  the 
tension,  T0  ,  as  a  function  of  the  pressure,  p  ,  the  reference 
radius,  rQ  ,  the  number  of  filaments,  n  ,  and  the  filament  angle, 
j30  ,  at  the  reference  radius.  From  Equation  (31)  ,  the  minimum 
radius  to  which  the  filaments  extend  is  given  by 


R" 


min 


=  X 


mm 


=  2  C' 


1 


where  the  reference  radius  is  chosen  such  that  aQ  =  0  .  Also, 
from  Equation  (32) ,  the  filament  becomes  tangent  to  the  meridian 
plane  for 


Xmin 


1 


2C2 


2  2 
Thus  for  C  <  1/2  ,  the  filament  forms  a  cusp,  while  for  C 

>  1/2  a  smooth  loop  is  formed.  For  the  special  case  (3Q  =  45° 
(C^  =  1/2) ,  the  complete  sphere  is  covered  with  filaments  over¬ 
lapping  severely  at  the  poles.  The  filament  curve  is  obtained 
from  Equation  (26) 


and  the  tension  in  the  filaments  by 
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Note  that  T  vanishes  for  R2  =  2C2  -  1  ,  which  corresponds  to 
the  minimum  radius  determined  above  when  C2  >  1/2  .  For  the 
three  distinct  cases  C2  <  1/2  ,  C2  =  1/2  ,  C2  >  1/2  ,  the  fila¬ 
ment  angle,  p  ,  and  the  tension,  T/T0  ,  are  plotted  versus  the 
radial  coordinate,  R  ,  in  Figure  5.  Corresponding  sketches  of 
the  net  are  also  shown  in  Figure  5. 

Example  2  -  Spinning  Cone  With  Uniform  External  Pressure.  - 
Take  p  to  be  the  uniform  surface  pressure  on  the  outside  of  the 
cone  (p  <  0)  .  Then  p/pQ  =  1  ,  a  =  a0  ,  and  K  <  0  . 

Hence  from  Equation  (22) 


Q  (X) 


seed 


The  boundary  condition  on  Q  is  obtained  immediately  from  Equa 
tion  (28) 


Q  =  C2  -  K  seed  J 1  -  C2 


The  special  case  of  radial  filaments  emanating  as  linear  rays 
from  the  apex  along  the  cone  is  given  by  C  =  0  .  Then 

Q  =  -  K  seed 

For  this  case,  each  filament  is  in  equilibrium  with  only  pressure 
and  inertia  loading,  the  interfilament  forces  vanishing.  For 
C  /  0  ,  the  maximum  radius  is  given  by  Equation  (31)  and  K  can 
be  eliminated  by  the  use  of  Equation  (28) 


r2  =  x  =  1-1  J 1  -  C2  cosd  =  1  +  HI - ' > 

max  max  k  "  -  C2 

The  condition  that  the  filaments  be  tangent  to  the  meridian  plane 
is  given  by  Equation  (32)  and  merely  reproduces  the  condition  for 
radial  filaments  given  above:  Q  =  -  K  seed  .  The  filament 
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geometry  and  tension  are  obtained  by  substituting  the  expression 
given  above  for  Q (X)  into  Equations  (26)  and  (27). 


sec  3  = 

ox  +  ro2x2 


+  4  (1-C2) -4  \  1-C2  KsecQ.  (1-2X) 


2  \  l-c2  -  Ksecoc  (1-X) 


OX  +  1_  [q2x2  +  4(i_c2)  -4  -^1-C2  Kseca  (1-2X)  +  K2sec2a (1-4X  +  3X  )  ] 
2  “  2 

These  results  are  shown  on  Figure  6  for  particular  values  of  C  , 
a  and  K  .  For  the  parameters  chosen,  T/T0  is  maximum  at 
X  =  Xmax  .  Because  |3  =  tt/2  at  X  =  Xmax  ,  there  is  no  edge 
load  at  the  aft  end  of  the  cone,  however,  there  is  a  concentrated 
load  at  the  apex,  D  ,  given  by 


D  =  n  T  coset 


where  n  is  the  total  number  of  filaments. 

Example  3  -  Spinning  Paraboloid.  -  Q  ^  0  ,  K  =  0  .  The 
meridional  curve  for  a  paraboloid  with  aperture-to-focal-length 
ratio  of  2  is 


Z  =  -  R2/4 


By  Equations  (22)  and  (28) ,  the  angular  velocity  parameter,  0  , 
is  related  to  the  boundary  condition  on  p  at  the  outer  edge  of 
the  surface 


0  =  1-  cos2  P0  cos2a0  =  4-^— 
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where  C  =  sinpQ  ,  and  the  value  of  cos2a0  =  1/5  is  obtained 
from  the  equation  of  the  meridian  and  simple  differential 
geometry.  The  filament  tension  and  geometry  in  the  surface  are 
easily  obtained  from  Equations  (27)  and  (26)  respectively. 


T  _ 

To 

1 1 
2  1 

x  4  +  C2 
5 

+  \x2  4  +  C2  +  4  (1  -  c2) 

■  5  5 

/  \  % 

cosp 

2 

^1  -  c2 

(-  +  i) 

\5X  5/ 

X 

(Hr1) 

+  Jx2/4_±_£i)+  i  (1  .  C2) 

These  results  are  shown  on  Figure  7  where  both  T/TQ  and  3 
are  plotted  versus  R  with  the  square  of  the  sine  of  the  filament 
angle  at  the  outer  edge,  j30  ,  as  a  parameter.  The  minimum  radius 
to  which  the  surface  is  covered  decreases  to  zero  as 
approaches  tt/2  .  At  this  limit,  however,  the  solution  is  a 
family  of  disconnected  parallel  circles. 

Except  for  the  limit  case  of  =  tt/2  ,  tensile  forces  exist 

at  both  the  inner  and  outer  edges  of  the  surface  *  These  forces 
are  obtained  from  the  relation  between  Q  and  (30  above  and  the 
definition  of  fi  in  Equation  (12) . 

n  =  m'^2ro2  =  4  +  C2 


Thus 


5  m'uu2r  2 
T°  4  +  sin2  (30 


This  is  the  tension  in  the  filaments  at  R  =  1  or  r  =  rQ  ,  the 
outer  radius  of  the  paraboloid.  The  edge  load  at  the  inner  edge 
is  obtained  from  Figure  7  and  the  value  of  TQ  above. 
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Solution  When  A  ^  0  .  -  For  the  general  loading  conditions 
A  ^  0  ,  0^0,  K  ^  0  ,  Equations  (19)  and  (20)  can  be  reduced 
to  one  algebraic  equation  and  one  differential  equation.  With 
minor  modification.  Equations  (19)  and  (20)  can  be  written  as 


d(T/T0  sinp)  2 
dR 


d  (T/T0  cos  (3)  2  =  2  T/T0  R  K  cos  (3  +  2  T/Tq  A  _  2  (T/Tq  cosp)2  d(cosg) 
dR  cosa  sina  cosoc  cosa  dR 


d (T/T0) 
dR 


+  2Q  T/Tq  R-2A  T/Tq  etna  -  1  (T/T0  sinP)2 

R 


If  these  expressions  are  added  together,  the  terms  involving 
derivatives  of-  T/T0  drop  out  and  the  resulting  expression  may  be 
solved  for  T/Tq  .  Thus, 


vlR.2  +  +  a  R  tana 

T  _  _ cosa _ 

T0  sin2p  +  cos2!  d (cosa) 


(33) 


Equation  (33)  in  conjunction  with  either  of  Equations  (19)  and 
(20)  or  some  combination  thereof  may  be  used  to  numerically 
determine  the  filament  tension,  T/T0  ,  and  the  filament  angle, 

P  ,  for  the  given  geometry  and  load  conditions. 

Example  4  -  Gravity  Loaded  Paraboloid.  -  A^0,  Q=K=0 

As  a  simple  example  solution  when  A  ^  0  ,  consider  a  paraboloid 
supported  at  its  outer  edge  in  a  gravitational  or  other  accelera¬ 
tion  field.  The  meridional  curve  is  again  as  in  Example  3 
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Equation  (33)  with  Q  =  K  =  0  reads 


AR  tana 


•  2  a  .  cos2  B 
sm^  B  +  - H 


cosa 


R 


d  (cosa) 
dR 


which  for  the  paraboloidal  geometry  may  be  written 


T  =  2A  (4  +  R2) 
TQ  4  +  R^  sin^B 


The  second  equation  used  is  obtained  by  eliminating  A  between 
Equations  (19)  and  (20)  and  then  carrying  out  the  indicated  dif- 
ferentations  and  solving  for  dB/dR  . 


d£ 

dR 


4RZ  cos 


(4  cscB  -  8  sinB  -  R^  sinB)  -  (4  +  R^)  ^  (4  +  R^  sin^B)  tanB 


4R  [R2  (1  +  cos2  B )  +  4]  (4  +  R2) 


A  digital  computer  program  was  written  and  used  to  determine 
T/T0  and  B  from  these  two  governing  equations.  The  initial 
condition.  Bo  at  R  =  1  ,  determines  the  minimum  radius  to  which 

the  solution  extends.  Thus  by  trial  and  error,  the  value  of  Bo 
for  which  the  paraboloid  is  completely  covered  may  be  determined. 
For  a  value  of  B0  ~  15.7172°  ,  the  solution  extends  to  R  < 

0.003  .  For  values  of  B0  <  15.7172°  ,  B  =  0  at  the  minimum 
radius  to  which  the  solution  extends,  thus,  a  force  is  required 
at  the  inner  edge.  For  Bo  >  15.7172°  ,  B  =  tt/2  at  the  inner 
edge,  thus,  the  filaments  end  only  at  the  outer  edge  and  no  force 
is  required  at  the  inner  edge. 

The  tension  variation,  T/(T0A)  ,  and  filament  angle  B  are 

plotted  versus  radial  position,  R  ,for  B0  =  15.7172°  on  Figure 
8.  Note  that  T/ (T0A)  is  maximum  at  the  outer  radius  R  =  1 
and  by  use  of  the  definition  of  A  (Equation  12)  and  Figure  8 
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1 


Tmax 
m '  ar 


o 


2.455 


This  equation  can  be  used  to  determine  the  maximum  size  paraboloid 
that  can  be  made  in  terms  of  the  filament  density  and  maximum 
stress.  In  engineering  terminology 

?max  _ 
p  n  2.455  °max 

where  p  is  the  weight  density,  n  is  the  number  of  g's  accel¬ 
eration  and  a/p  is  termed  the  specific  strength. 

The  Indirect  Problem  -  Isotensoid  Design 

Return  now  to  the  governing  Equations  (14)  through  (18)  and 
stipulate  that  the  filament  tension  be  uniform.  Then  T/Tq  =  1 
and  Equation  (14)  yields  an  explicit  expression  for  the  inter¬ 
filament  force 

^  =  Q  R  ctnp  sina  -  A  ctnp  cosa  (34) 

Equation  (34)  can  be  used  to  eliminate  Fc  from  Equation  (16)  and 
for  T/Tq  =  1  >  Equation  (15)  and  (16)  read 


d  (cos {3  cosa)  =  KR  p/p  +  A  secp  csca  (35) 

dR  ° 


d(R  sin_g_)  _  q  r2  cscq  _  AR  etna  csc|5  (36) 

dR 

A  convenient  form  for  these  equations  is  obtained  by  multiply¬ 
ing  Equation  (35)  by  cosg  cosa  and  Equation  (36)  by  R  sing  . 

Then 
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d(cos  p  cos2aj  _  2 (k  R  p/pQ  cosP  cosa  +  A  etna) 
dR 


(37) 


d(R2  s,in,2.§.)  =2  (fi  R3  -  A  R2  etna) 
dR 


(38) 


As  with  the  direct  problem,  the  presence  of  A  ^  0  compli¬ 
cates  the  integration  of  Equations  (37)  and  (38) .  As  a  conse¬ 
quence,  when  A  ^  0  ,  numerical  integration  has  been  employed  in 
their  solution.  Expressions  for  da/dR  and  dp/dR  suitable  for 
numerical  integration  are  readily  obtained  from  Equations  (37) 
and  (38)  by  carrying  out  the  indicated  differentiation.  There 
results 


da  _ 
dR 


-1 


-  K  p/pQR  csca  seep  -  Q  R etna  seczp  -  A  sec  p+  tan  p  ctnaR  (39) 


M  = 

dR 


Q  R  cscP  seep  -  A  etna  esep  seep  -  tanpR 


-1 


(40) 


and  from  Equation  (18) 


dZ 

dR 


etna 


It  is  clear  that  as  a  approaches  zero,  Equation  (18)  becomes 
unsuitable  for  numerical  integration.  This  is  remedied  by  chang¬ 
ing  independent  variables  from  R  to  Z  .  Then  Equations  (39) , 
(40) ,  and  (18)  read 


p/p  r  seca  seep+fi  R  sec 
dZ  *o 


P  +  A  tanp  sec 
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(42) 


=  -QR  tana  cscp  seep  +  A  cscp  seep  +  tanP  tanaR' 
dZ 


dR 

dZ 


tana 


(43) 


Integration  of  Governing  Equations 


Solutions  to  the  governing  equations  are  given  below  in  order 
of  increasing  complexity,  starting  with  cases  where  only  one  load 
parameter  is  nonzero  and  moving  on  to  all  load  parameters  nonzero. 
The  following  cases  are  considered: 


Case  1 
Case  2 
Case  3 
Case  4 
Case  5 


K  ?  0 
K  =  0 
K  =  0 
K  ^  0 
K  £  0 


n  =  0 
n  ^  o 

Q  =  0 
fi  >  o 
0  >  0 


A  =  0 
A  =  0 
A  £  0 
A  =  0 
A  ^  0 


Case  1  -  K  /  0.  -  Pressure  load.  This  case  is  treated  exten¬ 
sively  in  Reference  8. 

Case  2  -  fl  ^  0.  -  Load  due  to  spin.  For  this  case.  Equations 
(37)  and  (38)  read 


(a) 


d(cosP  cosa)  -  q 
dR 


(b) 


(UR2. . S^2J-)  =  2  a  R3 

dR 


Using  the  previous  notation  for  fiber  angle  at  the  reference  radius 
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sinP(R=l)  =  C 

where 

a (R=l)  =  0 

and  using  an  auxiliary  parameter 


expressions  (a)  and  (b)  can  be  readily  integrated,  yielding 


Substituting  (d)  into  (c)  yields  an  expression  for  the  meridian 
angle 


and,  using  Equation  (18) 


Expression  (f)  can  be  integrated  as 
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Other  quantities  of  interest  are  the  fiber  azimuth,  from  Equation 
(7)  and  from  expressions  (c)  and  (d)  above: 


/h)  df_  _  1_  sing  =  1  J  R4  +  (B-l) 

dR  R  cosg  sind  R  yj  (1_r2)  (r2-b+1) 

The  nondimensional  fiber  length  dL  =  d-t/ro  from  Equation  (2) 
and  (d)  becomes 


or,  with  Z  as  a  variable 


(j) 


dL  =  1 

dz  cosp  cosa 


=  cscp 


(R=D 


Expression  (c)  permits  the  calculation  of  fiber  lengths 
between  fiber  intersections,  defining  the  metric  properties  of  the 
two-family  fiber  net.  Expression  (j)  shows  the  fiber  length  incre¬ 
ment  in  the  axial  direction  to  be  dependent  only  on  the  parameter, 

C  ,  (i.e.  the  fiber  angle  at  the  equator  R  =  1  )  and  therefore 

relates  the  fiber  length  (and  consequently  the  structural  mass) 
directly  to  the  axial  distance  between  poles  of  the  spinning  struc¬ 


ture.  Hence,  with  the  boundary  condition 


2  (R=l) =  ° 


33 


An  examination  of  (g)  and  (d)  reveals  the  existence  of  real 
solutions  for  either  positive  or  negative  values  of  Q  .  Using 
the  boundary  conditions  Z  =  0  at  R  =  1  ,  (g)  becomes  for 

positive  values  of  Q 


(44a) 


And  for  negative  values  of  0 


Meridional  shapes  and  fiber  patterns  are  shown  in  Figures  9  and  10, 
respectively.  They  can  be  described  as  follows : 

1  For  0  <  C1 2  <-  Q  (or  0  <  B  <  1)  ,  the  structure  forms  a 

2 

barrel  open  at  both  ends,  symmetrical  to  the  equatorial  mid-plane 
at  R  =  1  with  the  fibers  terminating  at  the  barrel  ends  at  a 

radius  R  =  -Vr-TT  and  in  directions  tangential  to  the  meridian. 

2.  For  C2  =  (or  B=l)  ,  the  meridian  of  the  structure 
will  intersect  the  Z-axis  twice,  with  a  slope,  approached  in  the 
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,  and  the  axial  distance 


limit  by  cosa 


=  +VI^ 


,2  _ 


sin  (3 


(R=0)  '  r  (R=l) 

between  the  two  intersections  is  in  nondimensional  notation 


(2zr-o> 


The  structure  is  that  of  a  spindle,  symmetrical  with  respect 
to  an  equatorial  plane,  pointed  at  two  ends,  covered  by  fibers  that 
form  a  maximum  angle,  gQ  ,  at  the  equator  and  with  the  fiber  angle 
declining  to  zero  at  the  two  apexes.  This  structure  has  the  addi¬ 
tional  interesting  property  that  the  distances  between  fiber 
intersections  are  all  equal  and  given  by 


This  is  seen  readily  from  Expression  (i)  by  letting  B  =  1  and 
and  A$  =  2  .  The  fiber  pattern,  projected  into  the  equatorial 

plane  assumes  the  form  of  circles  with  radius  p  =  ^2.  tangent  to 

the  periphery,  R=1  ,  and  intersecting  the  center  (see  Fig.  11) . 

3.  For  ~  <  <  Q  (or  1  <  B  <  2  ) ,  the  structure 

describes  a  periodically  ncorrugatedn  tube  of  radius  range 
1  ^  R  >  ^ B-l  indefinitely  extended  along  the  Z-axis.  The  fiber 
angle,  p  ,  is  at  a  maximum  at  both  radial  extremes.  The  axial 
distance  between  two  extremes  is,  in  nondimensional  notation 


2z(-w  Vi  (i-c2» 

2 

4.  For  C  =  Q  ,  or  (B=2)  ,  Equation  (44)  yields  an 

indefinite  value,  real  only  for  R=1  .  The  corresponding  structure 
is  a  right  circular  cylinder  of  radius,  R  ,  covered  by  fibers 
forming  helices  of  constant  pitch. 


► 
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5.  For  C2  >  ft  or  (B  >  2)  ,  the  structure  is  geometrically 

similar  to  that  described  in  (3),  but  the  radius  range  is  restrict¬ 
ed  to  J(B-T)  >  R  >  1  .  It  can  be  readily  shown  that  the  two  types 
of  solutions  differ  only  by  a  geometrical  scale  factor  of  yj  (B-l) . 

6.  For  C=1  ,  a  special  case  obtains:  (h)  yields  Z  =  0 

for  r  <  1  The  structure  is  that  of  a  flat  annular  disk  origi¬ 
nally  discovered  by  Kyser  (Ref.  11) » 

7.  For  fi  <  0  (B  <  0  since  2C2  is  positive  definite) 
real  solutions  for  the  surface  exist  only  for  R  >  1^  andjreal 
solutions  for  the  filament  paths  exist  only  for  R<  \1-B  .  Thus 
for  0  <  0  (B<0)  ,  solutions  exist  in  the  range  1<RS  \1-B  .  The 
structure  is  formed  by  inward  arches  tangent  to  the  cylinder  R— 

At  the  outer  perimeter,  the  fibers  become  parallel  to  the  meridian, 
i.e.,  fiber  angle  {3  becomes  zero.  The  slope  angle  of  a  meridian 
in  the  surface  at  the  outer  periphery  is  equal  to  the  fiber  angle 


at  R=1 


l .  e . 


=  P  3. 

K  %ax 


=  3 


R=1 


Case  3-A^O  q  =  k  =  0 .  -  Uniform  acceleration  in  negative 
Z-direction.  For  this  case.  Equations  (37),  (38),  and  (18)  read: 


(1) 


d  (cos2  ft  cosja) _  2a  ctna 
dR 


d(R2  sin2!) =  _2A  r2  ctna 
dR 


(m) 


dZ 

dR 


ctna 


Eliminating  ctna  from  (1)  by  use  of  (m)  and  integrating 


cos2p  cos^a  —  -  2AZ  + 


(45) 


where  is  an 

of  Equation  (45) 


integration  constant.  Note  that  the  right  side 
is  confined  to  the  interval  (0,  1)  .  Equation 
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(45)  allows  the  determination  of  the  shape,  Z(R)  ,  once  a  and 
P  are  known.  These  functions  are  most  conveniently  determined 
numerically  using  Equations  (39)  and  (40)  which  for  0  =  K  =  0 
read 


(n) 


—  =  -  A  sec^(5  +  tan^p  etna  R 
dR 


(o) 


= 

dR 


-  A  etna  esep  seep  -  tanp  R 


-1 


The  shape  Z(R)  is  determined  by  Equation  (45).  These  Equations 
have  been  numerically  integrated  and  the  resulting  shapes  Z(R) 
are  shown  on  Figures  11  for  various  ranges  of  aQ  ,  p  ,  rQ  ,  and 
A  .  ° 


The  meridional  shapes  shown  in  Figures  11  terminate  at  both 
peripheral  edges  when  p  ,  the  filament  angle,  goes  to  zero.  The 
direction  and  magnitude  of  the  edge  loads  required  to  support 
these  surfaces  is  determined  by  the  surface  angle,  a  ,  at  the 
periphery,  the  number  of  filaments,  and  the  filament  load.  Recall¬ 
ing  the  definition  of  the  axial  load  parameter  (Eq.  12) ,  A  is 
the  product  of  the  reference  radius  and  the  ratio  of  surface  load, 
m'a  ,  to  filament  tension,  T-  .  Thus,  considering  rQ  constant 
as  A  approaches  zero,  the  filament  tension  dominates  the  surface 
load  and  the  shape  is  determined  almost  entirely  by  the  edge  loads. 
In  the  limit,  A  =  0  ,  there  is  no  surface  load,  and  the  surface 
becomes  a  hyperboloid  of  one  sheet  generated  by  the  now  straight 
filaments  (See  Flugge,  Ref.  3,  p.  78) .  This  is  indicated  in 
Figure  11a  where  and  pQ  are  maintained  constant  and  the 

meridional  shapes  shown  reflect  a  parametric  variation  of  A  . 

For  small  values  of  A  ,  Figure  lib  demonstrates  that  the  effect 
of  PQ  on  the  meridional  shapes  is  to  change  the  near-asymptotic 
angle  of  near-hyperbolic  curves.  Figure  11c  demonstrates  that 
the  isotensoid  generated  by  these  conditions  are  circularly 
symmetric  arches  supported  on  two  peripheral  lines  by  forces 
with  both  axial  and  vertical  components.  If,  however,  the 
upper  portion  of  the  surface  is  terminated  where  a  horizontal 
tangent  occurs  (a  =  tt/2)  the  axial  surface  of  load  A  is  sup¬ 
ported  entirely  at  the  base  periphery  and  the  upper  periphery 
requires  only  radial  support  loads.  A  wire  model  of  this  con¬ 
figuration  was  constructed  and  is  shown  in  Figure  12 .  The 
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intersections  of  the  wires  are  spot  welded  to  allow  interfilament 
force  transfer.  Note  also  that  a  simultaneous  reversal  of  the 
sense  of  the  acceleration,  "a"  ,  and  the  filament  load,  T0  ,  does 
not  affect  either  the  load  parameter,  "A"  ,  or  the  solution  shapes 
Thus  the  surfaces  described  by  the  meridional  curves  in  Figures  11 
are,  for  axial  loading  conditions,  either  isotensoid  or  isocompres- 
soid  surfaces  depending  on  the  surface  orientation  with  respect  to 
the  surface  load  and  the  edge  condition.  The  model  shown  in  Fig¬ 
ure  12  is  isotensoid  if  it  is  suspended  in  a  gravity  or  accelera¬ 
tion  field  from  the  smaller  internal  periphery  with  the  large  edge 
down,  or  an  isocompressoid  "cantilever  dish"  if  it  is  supported  at 
the  inner  edge,  with  the  inner  edge  down.  Application  of  this 
cantilever  dish  to  civil  engineering  structures  is  indicated. 

Figure  lid  shows  the  meridians  of  the  isotensoid  (or  isocom¬ 
pressoid)  shapes  generated  by  axial  loading  for  large  values  of  rQ  . 

An  interesting  and  practical  shape  is  obtained  as  the  radius 
r  approaches  infinity.  To  take  this  limit,  let  r  =  rro  +  y 
where  r^,  is  a  parameter.  Then  Equations  (n)  and  (o)  read 


(q) 


dB-  _  _  A  etna  _  tanB 
dy  sinB  cosB  r^  +  y 


and  Equation  (18)  becomes 


(r)  —  =  -  etna 

dy 

where 

A  =  A/rQ  =  m'a/T0 
A  can  also  be  written  as 
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Now  let  r^  -*  «  .  (y,  z)  are  now  cartesian  coordinates  and 

Equations  (p) ,  (q) ,  and  (r)  read 


da  _  p 
—  =  _  A  Sec  0 


A  etna 
sinp  cosj3 


etna 


This  system  of  equations  is  easily  integrated  with  the  follow¬ 
ing  result: 


sin  (Pmax) 

Z  =  — - : - — -  1  - 

2  A 


2  A 

sin  (2  R 


sina  cosp  =  cosB 


etna  = 


COS'1  R 


cos  {3  =  COSp 


;*( 


)T 
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The  integration  constants  were  evaluated  such  that  Z  0  at  the 
maximum  value  of  y  ,  and  by  Equation  (47)  j3  -  Pmax  a  ~  . 

The  values  of  maximum  z  and  y  are 


zmax 


y 

■^max 


sin^Pmax  = 
2A 


sin2  p 


max  _ 


2A 


sin2  p 


max 


2g 


sin2  p 


max 


2g 


(50) 


(51) 


Equation  (46)  defines  a  strikingly  simple  isotensoid  shape; 
a  cylindrical  parabola.  This  uniform  stress  shape  can  be  used  in 
civil  engineering  structures.  The  acceleration  is  then  the  gravi¬ 
tational  field  g  =  1  ,  and  the  internal  loads  are  uniform 
compression. 

The  surface  shapes  determined  from  Equation  (46)  are  shown  in 
Figure  13a.  Here,  zK  is  plotted  versus  yA  ,  thus  the  complete 
range  of  shapes  is  given  by  parametric  variation  of  Pmax  °nly* 
The  filament  angles  determined  from  Equation  (49)  are  shown  in 
Figure  13b.  Here  p  is  plotted  versus  yA  and  once  again  the 
complete  range  of  variations  is  given  by  parametric  variation  of 
p  .  Because  of  the  way  in  which  the  design  parameters  pmax 
and*  A  enter  Equation  (46) ,  the  ratio  of  zmax  to  ymax  and 
the  focal  length  can  be  varied  independently .  In  fact  from  Equa¬ 
tions  (50)  and  (51) , 


zmax 

^max 


tanp 


max 


2 


1 


and  from  standard  relations  between  focal  length  and  the  coeffi¬ 
cients  in  Equation  (46) , 


Focal  Length 


2d 
COS  P 


max 


2A 
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A  model  of  this  structure  was  made  with  slotted  aluminum 
strips  in  order  to  demonstrate  the  foldability  of  the  configuration. 
Pictures  of  the  model  in  hoth  the  deployed  and  folded  configura¬ 
tions  are  shown  in  Figures  14. 

Case  4  -  K  /  0,  Q  ^  0|  A  =  0.  -  In  the  remainder  of  this  re¬ 
port  the  applied  pressure  is  assumed  uniform,  p  =  PQ  •  With  this 
stipulation  and  A  =  0  ,  the  governing  Equations  (35)  and  (36) 
read : 


d(cosj3  cosa)  _  KR 
dR 


(52) 


d  (R  sing)  _  Q  R' 


dR 


(R  sing) 


(53) 


Equation  (52)  can  be  integrated  directly,  resulting  in 


cosa  =  [KR2/2  +  C,  ]  [cosg]  1 


(54) 


Equation  (53)  can  be  readily  integrated,  resulting  in 


R  sin|3  = 


(55) 


where  and  C2  are  constants  of  integration.  With  the  con¬ 

stant,  C2  ,  evaluated  such  that  (3  =  gQ  at  R  =  1  (and  with  the 
notation  singQ  =  C)  ,  and  with  C±  evaluated  such  that  a  =  0 
at  R  =  1  ,  and  letting  X  =  R  ,  Equations  (55)  and  (54)  are  com¬ 
pactly  expressed  as 


sing  -  [-  X  +  (C2  -  £)  -1  (56) 

L2  2  _x. 
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cosa 


+ 


|  (1  -  X) 

(c2  -  §)  - 

2  X 


(57) 


For  Q  =  0  ,  Equation  (56)  reduces  to  the  equation  for  geo¬ 
desic  paths  on  a  surface  of  revolution  and  Equation  (57)  reduces 
to  Equation  (30b)  of  Reference  8. 

With  the  intrinsic  filament  geometry,  (3  ,  and  the  intrinsic 
surface  shape,  a  ,  known  through  Equations  (56)  and  (57)  as  func¬ 
tions  of  the  radial  coordinate,  X  ,  the  more  practical  design 
functions,  Z  and  $  ,  can  be  determined  by  integrating  Equations 
(17)  and  (18) .  The  filament  length  can  also  be  determined  by 
integrating  the  first  of  Equations  (2) .  Changing  to  the  X-nota- 
tion,  these  equations  read 

dz  _  ^  etna 

^  2  yfx 

d$  _  _  tang 

dX  2X  sina 


dl  _  csca  secg 

dX  2 


(58) 


(59) 

(60) 


Solutions  of  Equations  (58)  and  (60)  in  terms  of  elliptic 
integrals  of  the  first  and  second  kind  are  given  below.  Equation 
(59)  which  determines  the  longitudinal  position  of  a  filament  as 
a  function  of  radial  position  leads  in  general  to  hyper-elliptic 
integrals  and  no  attempt  is  made  to  record  the  closed  form 
solution. 

In  order  to  integrate  Equations  (58)  and  (60)  it  is  convenient 
to  define  a  function,  y  ,  as 


y2  =  X  cos2  g 


.  2 
sin  a 


(61) 
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Use  of  Equations  (56)  and  (57)  to  eliminate  0  and  a  from  Equa¬ 
tion  (61)  yields  the  following  expression  for  y (X) 


(C2-|)  (62) 


Equations  (56) ,  (57)  ,  and  (61)  can  be  used  to  eliminate  0  and 

cl  from  Equations  (58)  ,  (59)  ,  and  (60)  with  the  following  result 


(63) 


di  _  _  tang 
dX  2X  sina 


2X  y 


(64) 


dl 

dX 


csca.  sec0  _  ^  1 _ 

2y 


(65) 


To 

a  cubic 


integrate  these  equations,  note  that  y  has  the  form  of 
in  X  ,  which  expressed  in  factored  form  is 


y 


2 


2 

(1  -  X)  (X  -  X,)  (X  -  X  ) 
4  ^  J 


(66) 


where 
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Investigation  of  these  roots  and  Equation  (66)  shows  that  for 
real  values  of  y  ,  four  types  of  solutions  are  possible  for 
0  <  <  1  and  the  range  of  X  including  X  =  1  , 


Type  A  X3  <  X2  <  X  <  1 

Type  B  X3  <  1  <  X  <  X2 

(69) 

Type  C  X  <  1  <  X3  <  X2 

Type  D  X2  ,  X3  complex  conjugates  X  <  1 

The  solutions  of  Equation  (63)  for  the  shape  Z (X)  and  Equa¬ 
tion  (65)  for  the  filament  length,  t (X)  ,  are  found  in  Reference 

12  for  the  four  types  listed  above. 
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Type  C: 


with 


and 


Type  D: 


(82) 
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Within  each  solution  type  are  a  multitude  of  isotensoid 
shapes.  For  the  direct  problem  where  K  ,  0  ,  and  C  are  speci¬ 
fied,  the  solution  type  can  be  identified  by  calculating  the  roots, 
X2  and  X3  ,  from  Equations  (67)  and  comparing  with  Inequalities 

(69).  With  the  solution  type  identified,  the  meridional  profile 
is  determined  by  the  appropriate  equation  for  Z(X)  (Equation  70, 
74,  78,  or  82)  with  the  aid  of  tabulated  elliptic  integrals. 

Figures  15  and  16  are  generated  from  Equations  (67)- (69)  and  indi¬ 
cate  the  solution  type  which  governs  for  various  values  of  the 
load  parameters,  K  ,  Q  ,  and  the  design  parameter,  C  .  For 
example,  in  Figure  15,  the  boundaries  between  solution  types  are 
defined  by  C  =  constant  curves  in  K,  fi  space;  and  in  Figure  16 
the  boundaries  between  solution  types  are  defined  by  0  =  constant 
curves  in  K  ,  C  space.  Note  that  solution  type-C  occurs  over 
such  a  limited  range  of  parameters  that  the  thickness  of  the  lines 
in  Figure  15  obscures  the  region,  and  in  Figure  16,  the  regions 
of  solution  type-C  is  barely  discernible. 

In  order  to  identify  and  discuss  the  isotensoid  meridional 
shapes  generated  by  the  solution  equations  over  the  entire  range 
of  the  parameters  K  ,  0  ,  and  C  ,  it  is  advantageous  to  resort 
to  a  semigraphical  technique. 

To  this  end,  consider  Equation  (57),  the  equation  for  cosa, 
and  note  that  the  denominator  in  the  expression  is  equal  to  cosP. 
The  expression  for  cosa  completely  defines  the  meridional  pro¬ 
files,  and  the  denominator  (cosP)  completely  defines  the  filament 
geometry  within  the  surfaces  generated  by  the  meridians.  Let  the 
numerator  be  called  N  and  the  denominator  D  .  Then 


cosa 


+ 


£ 


i 

X 


(88) 


The  sign  characterizes  the  symmetry  of  the  meridional  profiles 

with  respect  to  an  equatorial  plane  (for  convenience,  the  plane 
Z  =  0)  and  stems  from  the  boundary  condition,  a(R  =  1)  =  0. 
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The  meridian  profiles  depend  on  the  three  parameters  C  ,  0  , 
and  K  .  However,  0  and  K  dominate  the  meridian  profile  to  a 
greater  extent  than  C  .  The  semigraphical  method  of  identifying 
the  solution  meridional  profiles  is  facilitated  by  the  fact  that 
fi  and  C  appear  without  K  in  the  denominator  of  Equation  (88) , 
and  K  and  C  appear  without  0  in  the  numerator.  The  method 
is  straightforward  but  its  execution  is  tedious. 

For  the  moment,  let  C  be  a  fixed  constant.  Then  N  and 
D  individually  can  be  plotted  versus  X  on  the  same  graph  with 
H  as  a  parameter  in  the  plot  of  D  ,  and  K  as  a  parameter  in 
the  plot  of  N  .  The  ratio  of  the  values  of  the  separate  families 
on  the  graph  for  fixed  X  represents  the  cosa  at  the  value  of 
X  . 


A  great  deal  of  information  is  immediately  available  by 
examination  of  these  plots.  For  example,  when  D  >  1  ,  no  real 
solution  exists  fof  P  ,  thus  a  bound  on  the  solution  profile  is 
established.  When  N  >  D  ,  no  real  solution  exists  for  a,,  thus 
a  second  bounding  condition  is  apparent  from  observation  of  N 


and  D  .  Further, 

by 

setting  N  = 

D  ,  there 

results 

J  2 

K 

(1  -  X)  =  + 

J  o 

T  2  q\ 

1 

u 

1 

.H 

2 

V1  -  2X  " 

c  -2 

X 

(89) 

Equation  (89)  is  the  same  cubic  equation  encountered  in  the 
elliptic  integrals  for  Z  ,  and  its  roots  (1,  X 2  ,  X3)  are  given 

by  Equation  (67).  As  a  consequence,  by  observing  where  N  =  D 
on  the  plots  of  N  and  D  versus  X  ,  and  comparing  with  the 
inequalities  (69),  the  solution  types  can  be  easily  identified. 

Other  properties  of  the  solution  meridians  are  also  evident. 
When  N  =  0  ,  cosa  =  0  ;  thus  a  equals  tt/2  and  the  meridian 
locally  attains  a  maximum  or  minimum  value  in  Z  .  Further  if 
N  is  everywhere  nonzero,  the  meridional  profile  has  no  maxima 
or  minima  in  Z  .  When  D  =  1  ,  not  only  is  a  boundary  in  X 
reached,  but  P  =  0  at  this  boundary.  Also,  since  X  =  , 

X  >  0  is  a  requirement,  and  by  stipulation,  the  coordinates 
are  nondimensionalized  such  that  X  =  1  is  always  a  limiting 
value  for  the  meridional  curves.  (One  exception,  X  ^  1  at 
both  boundaries,  occurs  in  what  is  later  termed  special  case  1 
of  region  3 ) . 
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The  numerator,  N  ,  is  a  family  of  straight  lines  with  slope, 

-  J  2 

All  these  lines  pass  through  the  point  (X  =  1  ,  N  -  K1  -  C  ) 

which  is  also  an  intersection  point  for  D  .  This  is  * 
quence  of  the  boundary  condition  a  -  0  at  X  1  , 
cosa  —  1  and  N  —  D. 

Solution  Types  for  Pressure  and  Spin  Loading 

The  denominator,  D  ,  is  sketched  in  Figures  17a-i.  Examin¬ 
ation  of  Equation  (88)  reveals  that  the  D^cur^es  always  pass 

through  the  four  points  (X=1,D  =  +  ^(1-C) 


and  (X  =  -1 


-  ±  V 


D  =  +  Vl  +  C^).  Changing  the  value  of  C  shifts  these  points 
on  the^  X  =  +  1  lines  but  does  not  significantly  alter  the  char¬ 
acter  of  the  curves.  The  sketches  identify  three  naturally  sepa¬ 
rated  regions  for  values  of  the  parameter,  0  ,  and  the  two  inter¬ 
esting  boundaries  separating  these  regions.  The  sketches  show  how 
the  distinctive  curves  for  each  group  smoothly  pass  through  the 
limiting  boundary  cases  and  change  their  basic  character.  I 
examining  Figured,  note  that  e  is  a  smal!  positive  number 
introduced  only  to  show  that  the  parameter,  «  ,  is  close  to 
boundary  case.  The  regions  and  boundaries  are 


Region  1 

Q  < 

0 

Figures  17a, 

b 

Boundary  Case  1 

a  = 

0 

Figure  17c 

Region  2 

o  < 

„  2 

a  <  2c 

Figures  17d, 

e 

Boundary  Case  2 

n  - 

DO 

o 

to 

Figure  17f 

Region  3 

2C  2 

<  Q  <  00 

Figures  17g, 

h 

Typical  curves  for  the  three  regions  and  two  boundary  cases 
are  sketched  superimposed  on  Figure  18.  In  order  to  illustrat 
the  conditions  under  which  different  solution  types  are  encoun-^ 
tered,  a  limited  number  of  N-curves  are  also  shown.  The  in 
section  of  the  N-  and  D-curves  identifies  the  roots  (1,  X3  ,  X 2 ) 

of  the  cubic  Equation  (89) ,  and  by  comparison  to  Equation  (69) , 
the  solution  type  is  known.  The  region  of  interest  m  these 
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curves  is  in  the  vicinity  of  X  =  1  . 
evident. 


The  following  facts  are 


Region  1  0  <  0  (Figure  17a) .  -  Only  solution  types  A  and  B 
are  possible.  For  solution  type  A,  X2  (Equation  67)  is  always 

the  minimum  value  to  which  the  meridian  extends,  and  a  =  0  at 
that  limit.  The  upper  bound  is  always  X  =  1  at  which  a  =  0 
also.  For  solution  type  B,  X  =  1  is  always  the  minimum  value 
to  which  the  meridian  extends,  and  the  maximum  value  is  defined 
either  by  (Equation  67)  or  the  value  of  X  >  0  at  which 

D  =  1  .  The  upper  limit  will  depend  on  N  in  this  case. 

Boundary  Case  1  0=0  (Figure  17c) .  -  Only  solution  types 

A  and  B  are  possible.  For  solution  type  A,  the  bounds  are  the 
same  as  when  Q  <  0  ;  however,  for  solution  type  B,  the  maximum 
value  to  which  the  surface  extends  is  always  X  which  can  be 
quite  large  for  small  values  of  K  .  2 

Region  2 — 0  <  Q  <  2c2  (Figure  17d) .  -  only  solutions  of  type 
A  and  B  are  possible  in  this  region;  X2  is  always  the  minimum 

for  type  A,  and  X2  is  always  the  maximum  for  type  B. 

Boundary  Case  2  ^  —  2C  (Figure  I7f) .  —  Only  solution  types 

A  and  B  are  possible  here.  Note  that  only  under  the  conditions 
n  2 

^  =  2C  is  it  possible  for  the  meridians  to  close  in  the  axis 
of  symmetry.  This  is  illustrated  by  the  N— curve  labeled  A  on 
Figure  17f.  For  solution  type  A  ,  either  X2  or  zero  is  the 

minimum;  and  for  solution  type  B,  X2  is  always  the  maximum. 

2 

Region  3 — 2c  <  0  (Figure  17h) .  -  in  this  region,  all  four 
solution  types  are  possible  depending  on  N  .  This  region  is 
covered  in  detail  below. 


Meridional  Profiles  for  Pressure  and  Spin  Loading 

In  each  of  Figures  19a-g  there  are  two  families  of  sketches. 
In  the  upper  family  are  sketched  one  typical  D—  versus  X— curve 
and  a  number  of  N-  versus  X-curves.  The  lower  family  of  curves 
are  sketches  of  the  meridional  profiles  generated  by  the  N-lines 
above.  Note  that  these  meridional  profiles  are  topologically 
equivalent  to  the  physical  meridional  profiles  but  are  distorted 
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because  they  are  sketched  as  functions  of  X  -  R  ,  not  R  . 
Identifying  numbers  relate  the  meridians  generated  by  each  H 
line  shown .  Number  1  always  corresponds  to  large  negative  K 
K  increases  as  the  identifying  number  increases,  and  the 
highest  number  on  each  figure  corresponds  to  a  large  positive 
value  of  K  .  The  solution  type  associated  with  each  K-line 
is  also  recorded  in  the  upper  sketch. 


Boundary  case  1,  0=0,  Figure  19a.  -  The  N-line  marked 

corresponds  to  K  <  <  0  ,  passes  through  zero,  and  intersects 
the  D-curve  at  X  =  X2  >  1.  Thus  the  slope  of  curve  1  in  the 

lower  sketch  goes  through  zero  and  then  «  ,  forming  the  meridian 
of  closed  loops  shown.  As  K  increases,  the  loops  increase  in 
size  and  X ■  increases  as  depicted  by  curve  2.  In  the  limi  , 

as  K  —  0  curve  3  is  approached.  The  surface  shape  corresponding 
to  this  meridian  is  a  hyperboloid  of  one  sheet  generated  by^ 
straight  filaments  having  edge  load  but  no  surface  load  (K 
see  Flugge,  Ref.  3  p.  78).  The  slope  of  the  D-curves  at  X  =  1 

are 

/d  [D)\  _  1  C2  -  £ 

\dx/  X  =  1  2  Vi  -  c2 


and  the  slope  of  the  N—  curves  are 


1 


K 

2 


As  can  be  seen  on  Figure  19a,  when  the  slope  of  the  N-curve  at 
X  =  1  is  less  than  the  slope  of  the  D-curve,  solution  type  B 
results  with  X  >  1.  Thus  for 


0 


i  c2  -  n 
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N-curves  typified  by  number  4  generate  meridians  which,  when 
rotated  about  the  axis  of  symmetry,  generate  periodically  corru¬ 
gated  tubes  with  limits  on  X  of  1  <  X  <  •  When 


K 

2 


1  C2  -  fi 


the  corrugations  pass  to  the  limiting  case  of  a  cylinder  (curve  5) . 
Note  that  this  is  also  the  transition  from  solution  type  B  to 
solution  type  A.  N-curve  6  generates  corrugations  once  again; 
however,  the  limits  on  X  are  now  ^  X  <  1.  N-curve  7  is  a 

unique  boundary  between  corrugations  and  loops-  Note  that  for 

curve  7,  N-*-  0  and  N -►  D  at  the  same  time.  Thus  p-*-rr/2 

and  an  inward  cusp  is  formed.  N-curve  8  indicates  that  a  goes 
through  tt/2  and  then  on  to  tt  ,  forming  a  loop  inwards  as  shown. 
N-curve  9  generates  the  meridian  of  a  torus  which  is  the  boundary 
between  the  two  types  of  looped  meridians  shown  (8  and  10) .  There 
is  no  closed-form  equation  for  the  values  of  K  and  C  which 
generate  this  useful  shape.  However,  the  condition  for  this  merid¬ 
ian  is  that  Z  =  0  at  X  =  X^  .  This  is  solution  type  A.  Conse¬ 
quently,  the  torus  conditions  are  obtained  by  setting  Z  =  0  in 

Equation  (70)  and  X  =  X^  .  In  this  case,  ^  =  tt/2  (see  Equa¬ 

tion  72  for  X  =  X2)  and  the  incomplete  elliptic  integrals  become 

complete  elliptic  integrals.  The  condition  for  a  torus  is  there 
fore. 

(x  - 1  K  (ki>  - E  (ki»  - 0 

where  is  given  by  Equation  (73)  and  X2  and  X^  by  Equation 

(67)  with  Q  =  0  .  For  larger  values  of  K  ,  loops  of  type  10  in 
Figure  19a  result. 

2 

Region  2  0  <  fl  <  2C  ,  Figure  19b.  -  In  this  region  the 
D-curve  is  closed  and  D  is  always  less  than  1.  Curve  1  does 


53 


not  differ  from  the  case  fl  =  0  ;  however,  as  K  increases  from 
large  negative  values,  a  value  is  reached  where  Z  =  0  at 
X  =  X2  .  This  corresponds  to  curve  2,  a  toroidal  shape  with 

X  >  1  .  The  torus  is  the  boundary  between  the  two  types  of  loops, 
curves  1  and  3.  The  torus  conditions  are  obtained  by  setting 
Z  =  0  in  Equation  (74)  and  X  =  X2  .  As  K  increases  further, 

the  unique  case  labeled  "4"  occurs.  Here  D  and  N  both  go  to 
zero,  thus  an  outward  cusp  is  formed  with  3  —  tt/2  at  the  outer 
periphery.  The  condition  for  shape  4  is  that  N  =  0  at  X  =  X2  . 

For  higher  values  of  K  ,  corrugated  tubes  are  formed,  with  the 
limiting  case  again  a  cylinder,  6  ,  at  X  =  1  .  The  higher  N- 
curves,  7  through  11,  are  similar  to  those  discussed  under  bound¬ 
ary  Case  1  and  introduce  no  new  shapes. 

2 

Boundary  Case  2,  fl  =  2C  ,  Figure  19c.  -  N-lines  1-6 
correspond  to  solution  types  B  and  generate  meridians  already 
discussed:  the  two  types  of  closed  loops  reported  by  the  torus, 

an  outward  cusp  in  which  P  =  tt/2  at  the  periphery,  and  a 
corrugated  tube  degenerating  to  a  cylinder  when  the  slopes  of 
the  N-  and  D-curves  are  equal  at  X  =  1  .  For  larger  values  of 
K  solution  type  A  is  encountered,  and  corrugated  tubes  inboard 
of  X  =  1  are  generated  (curves  7).  The  unique  character  of 

2 

the  solutions  generated  when  0  =  2C  is  found  by  examining 
curves  8-12  which  are  the  only  solution  types  which  close  on 

2 

the  axis  of  symmetry*  When  ft  =  2C  ,  the  D-curve  is  a  parabola 
which  passes  through  the  three  points  (X  =  ,  D'=  0)  ,  (X  =  1  , 

D  =  ^1  -  c2)  ,  and  (X  =  0  ,  D  =  1)  N-curve  8  is  that  N-curve 
which  also  intersects  (X  =  0  ,  N  =  1) .  Under  these  conditions, 

X  =  X  =  0  .  When  these  conditions  obtain,  examination  of  the 
2  3  2 

solution  equations  for  solution  type  A  yields  =  1  ,  ^  = 

and  because  F  (tt/2  ,  1)  =  K  (1)  =  00  ,  the  meridional  shape  gener¬ 
ated  by  curve  8  is  asymptotic  to  the  X  =  0  line.  The  surface 
shape  is  therefore  a  closed  pointed  bottle  of  °°  extent  in  Z  . 
The  relationship  between  K  and  C  at  which  this  shape  is  gen¬ 
erated  is  obtained  by  setting  N  =  1  at  X  =  0  .  This  yields 


K  =  2  (  yi  -  C2  -  1) 
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Skipping  for  the  moment  to  N-curve  12,  note  that  n  =  -1  at 
X  =  0  .  It  is  easy  to  establish  that  Z  -*■  -oo  as  x  -*■  0  for 


curve  12,  and  that  K  =  2  (1  +  V-  C2)  defines  this  condition. 
For  values  of  K  between  these  two  limits,  the  meridians  close 
on  the  axis  of  symmetry  and  P  -►  0  as  X  -►  0  .  Thus  in  order 
for  the  solution  shapes  to  close  on  the  axis  of  symmetry, 


=  2C 


£ 


2  (  Kl  -  C  -  1)  <  K  <  2(1  + 


r- 


C2) 


and  P  =  0  at  X  =  0  . 


(90) 


Within  this  range,  a  distinct  value  of  K  generates  curve  10  for 
which  coscc  =  0  at  the  axis  of  symmetry.  This  occurs  when  N  =  0 
at  X  =  0  or 


K  = 


2 


C 


2 


thus  for  2  (  Vi  -C  -1)<K<2  shapes  typified  by 

curve  9  are  generated,  and  for  Z  (  ijl  -  Q2  <  K  <  2  (1  +  -  c2) 

shapes  typified  by  curve  11  are  generated.  These  surface  shapes 
are  closed  bottles  requiring  supporting  structure  at  the  apices 
while  shape  10  is  self  supporting.  Values  of  K  above 

2  (1  +  ^/l  -  C  )  generate  the  close-loop  meridians  represented 
by  curve  13.  It  is  not  likely  that  closed  loops  of  the  other 
type  (Curve  9  on  Figure  19b)  or  tori  are  generated  for  this 
value  of  Q  in  combination  with  positive  values  of  K  . 

2 

Region  3  Q  >  2C  ,  Figures  19d,  19e#  19f.  -  Region  3  is 
complicated  by  the  fact  that  the  D-curve  has  an  inflection  point. 
The  position  of  this  point  relative  to  X  =  1  has  a  distinctive 
effect  on  the  solution  types  generated. 

2 

It  can  be  shown  chat  for  1  >  C  >  1/3  the  inflection  point 
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(— —  d  =  0 )  occurs  in  the  range  X  <  1  .  For  C  <  1/3  ,  it 

ax2  ' 

occurs  in  the  range  X  >  1  provided 


i.  vr : 


+  3c  <  n-  < 


-  4C  +  3C 


When  Q  >  2  ,  the  inflection  point  occurs  at  values  at  X  where 
D  is  greater  than  1  (X  <  X  minimum) .  Thus  three  distinct  cases 
arise,  hereafter  referred  to  as  special  cases.  In  special  case 
1  the  inflection  occurs  at  X  >  1  ,  Figure  19d.  The  conditions 

for  this  case  are  C  <  1/3  and  0  bounded  by  Equation  (91) . 
These  are  the  only  conditions  under  which  solution  type  C  is 
encountered  where  X2  and  X3  are  greater  than  1  >  X  (N-curve 

7C  in  Figure  19d) .  In  the  range  of  0  for  which  Figure  19d 
applies,  N-curves  1-5  are  similar  to  those  discussed  before. 
N-curve  7C  typifies  the  only  range  of  solution  type-C. 

Note  that  a  fifth  solution  type  is  possible.  It  occurs  only 
under  the  same  conditions  as  solution  type-C  and  generates  cor¬ 
rugated  tubes  (curve  7).  For  this  solution  type  1  <  X3  £  X  <  X2 

This  is  the  exception  identified  earlier  where  the  solution  is 
not  bounded  by  X  =  1  .  The  solution  equations  for  this  rare 
case  are  not  recorded  here  but  can  be  obtained  from  Reference  12. 

Of  particular  interest  is  the  fact  that  in  this  range 

O 

of  0  and  C  ,  the  D-curve  is  concave  upward  at  X  =  1  .  As 
a  consequence,  when  the  slope  of  the  N-curves  is  equal  to  the 
slope  of  the  D-curves  at  X  =  1  ,  solutions  for  X  >  1  and 
X  <  1  occur  for  the  same  values  of  all  parameters.  Thus  a 
corrugated  shape  can  fare  into  an  inward  arch  with  compatible 
filament  patterns  at  X  =  1  •  This  is  shown  by  the  extension 
to  curve  6  in  the  lower  sketch  of  Figure  19d.  When  the  inflection 
point  is  at  X  =  1  ,  which  is  the  condition  which  separates  spe¬ 
cial  cases  1  and  2,  the  corrugated  tubes  degenerate  to  a  cylinder 
when  D  and  N  have  the  same  slope;  thus  the  interesting  merid¬ 
ian  which  generates  a  cylinder  with  an  end  cap  identified  by  the 
numbers  6—7  on  Figure  19d  is  obtained. 
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Once  the  slope  of  N  is  greater  than  the  slope  of  D  at 
X  =  1  ,  the  lower  limit  of  the  solution  shapes  is  reached  at  the 
value  of  X  where  D  =  1  and  0  =  0  .  Thus  curves  6-11 
generate  meridional  profiles  with  central  openings  requiring 
support  of  the  filaments.  That  unique  curve  corresponding  to 

cosa  =  0  (a  =  at  the  same  value  that  D  =  1  (curve  9) 

requires  only  hoop  constraint  while  all  others  require  tension 
or  compression  support. 

The  limiting  value  of  K  for  these  solutions  is  obtained 
by  setting  N  =  -1  at  X  =  X2  (curve  11) .  For  values  of  K 

larger  than  this,  the  closed  loops  like  curve  12  are  obtained. 

It  is  doubtful  that  tori  and  the  other  type  of  loops  exist  for 
large  values  of  K  and  this  range  of  ft  .  With  reference  to 
Figure  17h,  it  is  apparent  that  solution  type  D  is  encountered 
whenever  the  N-lines  intersect  neither  the  D-curve  for  nega¬ 
tive  X  nor  the  D-curve  for  positive  X  except  at  X  =  1  . 

This  occurs  for  two  ranges  of  K  ,  one  negative  and  one  positive. 

In  a  second  special  case,  the  inflection  occurs  between  the 
value  of  X  where  D  =  1  and  X  =  1  .  This  is  shown  on  Figure 

2  2 

19e.  The  conditions  for  this  are  2C  <  ft  <  2  and  1  >  C  >  1/3 

or  2C2  <  ft  <  2  ,  C2  <  1/3  and  Equation  (91)  not  satisfied.  The 
curves  associated  with  all  numbers  but  8  are  similar  to  cases 
already  covered.  For  N— curve  8,  the  N— line  and  D— curve  are 
tangent  at  X  =  X2  =  X3  .  This  represents  a  branch  or  bifur¬ 
cation  point  from  which  two  isotensoid  meridians  are  possible  as 
shown  in  the  lower  sketch  8  in  Figure  19e. 

The  third  special  case  is  obtained  when  ft  >  2  and  is  shown 
in  Figure  19f.  There  the  inflection  point  in  D  falls  outside 
the  range  of  real  solutions.  All  the  profiles  except  those  num¬ 
bered  8  and  12  have  been  discussed  above.  In  cases  8  and  12, 

D  =  +  1  at  X  =  X2  and  N  =  1.  Thus  for  case  8,  a  corrugated 

tube  is  formed  in  which  0  =  0  at  the  inner  radius,  and  for 
case  12  loops  are  formed  for  which  0=0  at  the  inner  or 
minimum  radius. 

In  general.  Region  3  includes  solutions  of  types  A,  B,  and 
D  and  for  certain  special  conditions,  type  C.  The  minimum  radial 
coordinate  to  which  the  solutions  extend  is  bounded  by  the  value 
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of  X  where  D  =  1  over  most  of  the  range  of  K  applicable  to 
solution  types  A,  D  and  C. 

In  Figure  19g  the  isocompressoid  meridional  profiles  are 
shown  for  0  <  0  .  Recalling  the  definition  of  Q  (Equation  12) 
the  filament  load,  Tq  ,  is  necessarily  less  than  zero  if  0  <  0. 

For  T  <  0  ,  the  sense  of  p  is  reversed  in  the  sign  conven- 
o  o  3 

tion  for  K  .  The  N-lines  1,  and  7  through  13  generate  merid¬ 
ional  profiles  already  discussed,  in  each  case,  however,  the 
pressure  loading  is  reversed  and  the  internal  loads  are  com¬ 
pressive.  For  N-curves  2  and  6,  N  =  +  1  respectively  and 
D  =  1  .  Thus  curve  2  generates  outward  loops  tangent  to  the 
value  of  X  at  which  X  =  and  D  =  1  simultaneously.  In 

this  case  0  =  0  at  the  outer  periphery.  For  curve  6  a  corru¬ 
gated  tube  is  formed  tangent  to  the  same  X  value  at  its  outer 
periphery  and  with  £  =  0  .  Curves  3,  4,  and  5  are  inward  arches 
with  positive,  zero,  and  negative  slope  at  the  outer  periphery 
defined  by  D  =  1  .  The  special  case  when  K  =  0  (not  shown  on 
Figure  19g)  was  discussed  under  the  general  case  2  for  which 
K  =  A  =  OH^O  . 


In  order  to  illustrate  briefly  the  design  potential,  four 
spinning  pressurized  isotensoid  designs  are  given.  Figure  20 
shows  a  closed  spinning  isotensoid  pressure  vessel.  This  design 
closes  on  the  axis  of  symmetry  and  is  obtained  by  choosing 
2 

Q  =  2C  and  determining  C  by  setting  N  =  0  at  X  =  0.  This 
corresponds  to  N-curve  10  in  Figure  19c.  Figures  21  and  22 

2 

are  spinning  isotensoid  "tires"  with,  0  =  0°  (0  >  2c  )  and 
2 

0  =  90°  (Q  <  2C  )  at  the  inner  periphery,  respectively.  They 
correspond  to  N-line  8  on  Figure  19b  for  zero  0  at  the  hub, 
and  N-line  9  on  Figure  19d  for  0  =  90°  at  the  hub. 

Figure  23  depicts  a  design  of  a  spinning  net  with  internal 
(positive)  pressure  on  the  upper  portion  and  external  (negative) 
pressure  on  the  lower  portion,  representative  of  the  loading  on 
a  spinning  aerodynamic  decelerator. 

Case  5-K^0,Q^0,A^0.  -  An  analytic  solution  was 
not  found  for  this  general  loading  case  where  all  loads  are  non¬ 
zero.  The  equations  suitable  for  integration  by  digital  computer 
were  used.  These  are  Equations  (18),  (39),  (40)  when  a  is 

sufficiently  far  from  zero;  and  Equations  (41),  (42),  and  (43) 
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when  a  is  sufficiently  far  from  tt/2  .  There  are  five  signifi¬ 
cant  parameters  in  these  equations:  the  three  load  parameters 
K  ,  0  ,  and  A  ,  and  two  boundary  parameters,  3q  and  aQ  at 

r  =  r  .  (The  third  boundary  parameter,  Zq  ,  is  unessential 

because  it  serves  only  to  shift  the  surface  with  respect  to  the 
Z-axis.)  This  results  in  a  five-dimensional  space  to  explore 
for  significant  shapes  using  numerical  techniques.  As  a  conse¬ 
quence,  no  attempt  to  classify  shapes  was  made.  Instead,  a 
modest  attempt  was  made  to  determine  the  effect  of  a  uniform 
acceleration  load,  A  ,  on  the  shapes  determined  in  the  previous 
case  treated  where  only  K-  and  Q-loads  were  present.  The 
details  of  the  calculation  as  well  as  the  computer  program  used 
are  given  in  the  Appendix. 

In  order  to  obtain  a  reasonably  broad  coverage,  3  was 

2  1 

set  equal  to  45°  and  held  constant  so  that  C  =  ~  .  Shapes  were 

2  2  2  ^ 

determined  for  Q  <  2C  ,  Q  =  2C  ,  and  Q  >  2C  for  various 
values  of  K  and  for  A  =  0  .  For  each  one  of  these  shapes 
determined,  a  second  shape  was  determined  for  A  =  0.25  .  The 
plots  of  the  resulting  meridional  profiles  shown  in  Figures  24 
illustrate  the  sensitivity  of  isotensoid  shapes  to  axial  loading 
A  ^  0  . 


In  particular,  note  that  none  of  the  isotensoids  with  the 
A-load  present  extend  to  the  axis  of  symmetry.  This  leads  to 
examination  of  the  governing  equations  for  values  of  R  near 
zero  to  determine  if  there  are  any  isotensoid  surfaces  with 
A  ^  0  which  close  on  the  axis  of  symmetry.  The  details  of 
this  study  are  included  in  the  Appendix.  The  results  can  be 
summarized  as  follows: 

In  the  vicinity  of  R  =  0  ,  the  filament  angle,  3  ,  and 
the  meridional  angle  ,  a  ,  are  governed  by  the  approximate 
relations 


p'  = 


(92) 
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(93) 


a 


It  is  apparent  from  Equation  (92)  that  U  >  A  is  a  nec¬ 
essary  condition  for  closed  shapes. 

In  order  to  illustrate  surfaces  of  this  type.  Equations  (92) 
and  (63)  are  used  to  generate  values  of  a  and  P  for  small  r 
which  can  be  used  as  initial  conditions  for  the  governing  Equa¬ 
tions  away  from  the  singularity  at  R  =  0  .  The  results  are 
shown  in  Figure  25  for  A  =  1  ,  K  =  1  ,  and  a  range  of  A  ;  and 
in  Figure  26  for  0  =  1  ,  K  =  0  for  a  range  of  A  .  The  varia¬ 
tion  of  P  with  R  is  shown  in  Figure  27  for  one  case,  0  =  1  , 
K  =  0  ,  A  =  0.1  .  This  is  representative  for  all  cases  shown. 
Note  that  P  =  0  at  the  lower  opening  in  the  surface,  thus,  a 
tension  ring  (or  the  load  equivalent)  and  an  axial  force  must  be 
supplied  for  equilibrium.  There  is  an  indication  therefore  that 
this  analysis  may  be  applied  to  approximate  conditions  for  design 
of  a  high-altitude  (spinning)  parachute  of  the  type  discussed  in 
Reference  13.  The  edge  condition  is  then  supplied  by  payload 
weight  which  is  suspended  from  the  filaments. 
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CLOSURE 


The  results  reported  here  represent  a  step  toward  reducing 
net  structure  design  to  standard  methods.  Consideration  has  been 
restricted  to  axisyrmnetric  nets  subject  to  axisymmetr ic  but,  in 
general,  non-normal  surface  loads.  The  filament  geometry  is 
restricted  to  two  families  of  filaments  symmetrically  disposed 
within  the  surface. 

In  the  problems  where  the  exact  surface  shape  is  prescribed, 
the  tension  variation,  filament  geometry,  and  interfilament  force 
level  are  given  in  terms  of  load  levels  and  boundary  conditions. 
Also  the  region  of  the  prescribed  shape  over  which  a  net  struc¬ 
ture  can  be  used  to  react  the  loads  is  given  as  a  function  of  the 
design  parameter,  C  =  sinf3Q  . 

The  exact  shape  cannot  be  prescribed  if  isotensoid  design  is 
desired  in  two-family  systems;  however,  some  freedom  in  design  is 
possible.  When  the  load  conditions  are  pressure  and  spin  (K  ^  0 
Q  ^  0  ,  A  =  0)  the  designer  may  refer  to  the  classification  of 
shapes  in  Table  I  and  determine  values  of  the  design  parameter, 

C  =  sinP  ,  which  generate  isotensoid  shapes  best  suited  for  the 

particular  design  application.  With  K  ,  ft  ,  and  c  known,  the 
surface  shape  may  be  determined  from  the  appropriate  closed-form 
solutions  which  require  use  of  elliptic  integral  tables,  or  the 
digital  computer  program  given  in  Appendix  A  may  be  employed. 

When  the  load  condition  includes  A  ^  0  ,  use  of  the  digital  com¬ 
puter  program  is  mandatory.  If,  in  this  case,  it  is  desired  that 
the  surface  close  on  the  axis  of  symmetry,  the  method  for  gener¬ 
ating  initial  conditions  near  the  apex  as  outlined  in  Appendix  A, 
may  be  employed. 

It  is  clear  from  the  restrictions  placed  on  the  solutions 
obtained  that  a  number  of  problems  cannot  as  yet  be  handled.  The 
most  glaring  omission  for  symmetric  loading  is  the  ability  to  pre¬ 
scribe  structure  shape  and  obtain  isotensoid  design  solutions. 

Two  directions  can  be  taken  in  attacking  this  problem,  both 
requiring  an  increase  in  the  number  of  layers  of  filaments  in 
the  net.  First,  a  finite  number  of  layers  as  treated  in  Refer¬ 
ence  8,  with  the  addition  of  interfilament  force  transfer  between 
families  in  each  layer,  Fc  ,  and  interf ilament  force  transfer 
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between  layers,  ,  will  allow  isotensoid  design  when  shape  is 

prescribed.  Second,  passage  to  a  continuum  of  filament  layers, 
as  done  in  Reference  10,  may  be  attempted.  Non-uniform  pressure 
and  meridional  surface  loads  must  be  included,  however,  and  the 
non-geodesic  nature  of  the  filament  paths  must  be  accounted  for 
so  that  the  internal  meridional  and  circumferential  loads  are 
accommodated.  This  leads  to  a  pair  of  integral  eguations  for 
the  distribution  of  filament  density  as  defined  in  Reference  10. 

As  discussed  earlier,  physical  interpretation  of  the  results  of 
this  method  requires  approximation  by  a  finite  number  of  layers. 
It,  therefore,  appears  that  the  assumption  of  a  finite  number  of 
layers  at  the  onset  is  the  best  approach.  If  an  exact  solution 
covering  the  required  region  of  the  prescribed  shape  is  not  possi¬ 
ble,  the  difference  between  the  resulting  shapes  and  the  pre¬ 
scribed  one  may  be  minimized  with  respect  to  the  boundary  values 
of  the  angles  in  each  layer. 

A  second  serious  omission  is  the  ability  to  design  isotensoid 
filament  surfaces  in  the  presence  of  asymmetric  loading.  It  is 
apparent  that  isotensoid  design  is  not  possible  if  the  load  dis¬ 
tribution  varies  with  time.  However,  for  fixed  asymmetric  loads, 
say  of  the  type  which  lead  to  resultant  moment,  shear  and  torque 
at  the  surface  edges,  the  conditions  under  which  isotensoid 
design  is  possible  can  be  determined.  This  could  lead  to  signi¬ 
ficant  weight  savings  for  some  design  applications. 

Finally,  the  role  of  interfilament  force  level  has  not  been 
exploited  as  a  full  design  parameter.  The  weight  of  a  net  will 
be  affected  by  the  level  of  the  interfilament  forces  because 
matrix  material  or  knots  are  present.  Thus,  it  may  be  desirable 
to  use  the  maximum  value  of  interfilament  force  as  a  design  para¬ 
meter  in  much  the  same  way  as  C  =  sinPQ  is  employed.  Methods 

should  then  be  developed  to  have  the  interf ilament  force  as  a 
design  input.  Also,  isotensoid  shapes  could  be  classified  in 
terms  of  this  parameter  instead  of  C  =  sin|3 

o 


62 


V. 


Reference 
Meridional  Plane 


r 


Meridional  Plane 


Figure  1. 


Coordinates  and  Surface  Load 


2c  Fiber  Layer  with 

Complimentary  Sets 
at  P  and  -P 


2d  Single  Fiber  Be¬ 
longing  to  Both 
Sets  in  a  Layer 


Definition  of  Filament  Sets 


Figures  2. 


and  Filament  Layers 


ro  1=3 


1  2  R 


Figure  9.  Meridians  of  Spinning  Isotensoids 

(Isocompressoids)  $  =  45°  B  =  1/Q 
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Figure  10.  Fiber  Path  of  Spinning  Isotensoids 

(Isocompressoids)  =  45°  B  =  1/^ 
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Figure  11a.  Computed  Meridional  Shapes  for  Axial 
Load  (K  =  ft  =  0,  0) 
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Figure  11c.  Computed  Meridional  Shapes  for  Axial 
Load  (K  =  Q  =  0,  A  M) 
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Figure  12. 


Cantilever  Dish,  " Isocompressoid"  Structure 
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Figure  13a.  Parabolic  Shape  of  Isotensoid 
(Isocompressoid)  Cylinders 
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Figure  14.  Model  of  Isotensoid  (I socompressoid)  Parabolic  Cylinder 


o 

Figure  15.  Spinning  Isotensoid  Network  Structures 


Ranges  of  Solution  Types  (With  3  as 
Parameters)  ° 
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for  Different  Regions  of 


Figure  17.  Typical  Sketches  of 
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Figure  18.  Composite  Sketch  of  D  for  All  Regions  of  fi 
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Figure  19a.  Isotensoid  Meridional  Profiles  for  Boundary  Case  1 

A  =  0  K^O  0=0 
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Figure  19c.  Isotensoid  Meridional  Profiles  for  Boundary  Case  2 

A  =  0  K^O  n=  2 c2 


87 


89 


K=l;  C= 


£T. 

2  ' 


r  . 
mm 


=  0; 


(sing) r  =1; 
mi  “ 


n=  1.5 


20. 


Spinning  Isotensoid  Pressure  Vessel 


K=l. 125; 


.=  jT. 

2' 


min/r 


=1/3;  (sing)  =0;  «=1.519 


Figure 


Spinning  Isotensoid  "Tire"  Fibers 
Radial  At  Hub 
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Figure  22.  Spinning  Isotensoid  "Tire"  Fibers 

Tangential  To  Hub 
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Figure  23.  Spinning  Isotensoid  Net  Loaded  with  Discontinuous 

Pressure  Profile 
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RADIUS 


RADIUS 


Figures  24.  Computed  Meridional  Shapes  for  Cases  4  and 

5.  a  =  0°,  g  =  45°,  Load  Parameters  as 
o  o  2 

Indicated  ( 0  >  2C  ) 


K=4. 83 


A-0 


Z 


fi=2 

K=9 


A=0 


Computed  Meridional  Shapes  for  Cases  4  and 
5  a  =  0°f  3  =  45°,  Load  Parameters  as 

o  '  o  2  ' 

Indicated  (2  >  2C  ) 


Figures  24. 


Figure  25.  Computed  isotensoid  Shapes,  Closed  at 
Center.  Loads  as  Indicated 
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RADIUS 

Figure  26.  Computed  Isotensoid  Shapes  With  Spin  and  Axial  Load 
Closed  at  Center 


Beta  (Deg) 


APPENDIX 


NUMERICAL  INTEGRATION  OF  GOVERNING  EQUATIONS 


The  most  convenient  form  of  the  governing  equations  for 
numerical  integration  are  the  equation  sets  (39),  (40),  and  (18), 

or  (41),  (42),  and  (44).  They  read  (for  uniform  pressure): 


2  2  2 

^  =  -K  R  csca  seep  -OR  etna  sec  P  -  A  sec  P  +  tan  P  etna  R 
dR 


—  =  Q  R  cscP  secP  -  A  etna  cscP  seep  -  tanP  R 
dR 


HA-1) 


dZ 

dR 


-etna 


—  =  k  R  seca  secP  +  A  R 
dZ 


dp 

dZ 


dR 

dZ 


-Q  R  tana  cscP  secP 


-tana 


o  2  2 

sec  P  +  A  tanP  sec  P  -  tan  P  R 

+  A  cscP  secP  +  tanP  tanct  R 


HA-2) 


A  digital  computer  program  was  written  for  use  in  integrating 
equation  sets  (A-l)  and  (A-2)  using  a  second  order  Runge-Kutta 
method.  The  program  listing  is  given  at  the  end  of  this  appendix. 
In  using  this  program  to  determine  the  surface  profile  Z(R)  , 
the  load  parameters  K  ,  fi  ,  A  and  initial  values  for  a  ,  P  , 
and  R  and  Z  are  selected.  The  program  can  be  then  applied  for 
either  increasing  or  decreasing  increments  of  Z  and  R  •  “ 

a  approaches  zero,  equations  (A-2)  must  be  used  because  of  the 
etna  term  in  (A-l),  while  if  a  approaches  rr/2  Equations 

(A-l)  must  be  used  because  of  the  tana  term  in  (A  2) .In 
transferring  from  one  equation  set  to  another  in  the  middle  of 
a  run  the  last  print  out  from  the  computations  using  one  set 
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of  equations  is  used  for  initial  conditions  in  the  other  set  of 
equations.  in  this  manner,  surfaces  which  go  through  0,  tt/2 
and  cyclic  increases  of  these  angles  can  be  determined. 

For  surfaces  which  approach  R  =  0  ,  neither  set  of  equa¬ 
tions  is  satisfactory  for  digital  computation  near  R  =  0 

because  of  the  appearance  of  R  in  both  equation  sets.  In 
order  to  remedy  this  situation.  Equations  (A-l)  were  examined  for 
limiting  restrictions  as  R  0  .  These  were  chosen  in  prefer¬ 
ence  to  Equations  (A-2)  because,  as  will  be  shown,  a  must 
approach  tt/2  and  dR/dZ  therefore  becomes  unbounded. 

Limit  when  R  0  .  -  From  the  first  two  of  Equations  (A-l) 

it  is  apparent  that  tan0  must  approach  zero  with  at  least  the 
first  power  of  R  .  For  the  special  case  A  =  0  ,  no  further 
restrictions  are  placed  on  any  of  the  variables  because  singulari¬ 
ties  due  to  0  -►  0  linearly  with  R  are  always  cancelled  by  the 
appearance  of  R  in  the  numerator.  However  if  a  /  0  ,  a 
second  singularity  appears  in  the  second  of  equations  (A-l)  which 
can  only  be  cancelled  if  a -►  tt/2  as  0  0  .  Thus,  all  isoten- 

soid  shapes  have  zero  slope  (dZ/dR  =0)  at  R  =  0  for  A  /  0  . 
In  order  to  compute  these  shapes,  it  is  necessary  to  determine 
proper  initial  conditions.  To  this  end  assume 


13  -  C-L  R  +  C2  R2  R~0  (A-3 ) 

and 

a  =  rr/2  -  C3  R  -  c4  R2  R  -  0  (A-4) 


The  constants  are  evaluated  by  substituting  for  0 

and  a  into  Equations  (A-l) ,  making  appropriate  approximations 
for  small  R  and  matching  coefficients  of  powers  of  R  .  After 
some  manipulation,  there  results 


0  = 


(A- 5) 
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(A—  6) 


a  =  tt/2  -  AR 


K  r2 
2  R 


2  3 

AR  _  KR_ 

2  6 


(A-7) 


Equations  (A-5) ,  (A-6) ,  and  A-7)  can  be  used  to  generate  initial 

conditions  for  a  ,  p  ,  and  Z  which  generate  isotensoid  shapes 
closing  on  the  axis  of  symmetry.  Note  that  real  solutions  are 
only  possible  if 


n 


>  0 


(A-8) 


Aside  from  the  restriction  of  Equation  (A-8) ,  one  can  for 

2 

all  combinations  of  K  ,  0  ,  A  (Q  -  A  >0)  determine  initial 
conditions  which  generate  surfaces  which  close  on  the  axis  of 
symmetry. 

Illustrative  examples  of  Q  =  1  ,  K  =  1  ,  and  Q  =  1  ,  K  =  0 
both  with  parametric  variation  of  A  are  demonstrated  in  Figures 
24  and  25.  Figure  26  shows  the  fiber  angle  for  a  typical  loading; 
n=l,K=0,A=0.1. 
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